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ABSTRACT 



mixing angles can be generated in the context of a natural theory for the TeV-scale. Sneutrino- 
mixing phenomena provide valuable tools for connecting the physics of neutrinos and supersym- 
^ , metry. We examine the theoretical structure of the seesaw-extended MSSM, retaining the full 

' complexity of three generations of neutrinos and sneutrinos. In this general framework, new 

flavor-changing and CP-violating sneutrino processes are allowed, and are parameterized in terms 
of two 3x3 matrices that respectively preserve and violate lepton number. The elements of these 
matrices can be bounded by analyzing the rate for rare flavor-changing decays of charged leptons 
and the one-loop contribution to neutrino masses. In the former case, new contributions arise 
in the seesaw extended model which are not present in the ordinary MSSM. In the latter case, 
sneutrino-antisneutrino mixing generates the leading correction at one-loop to neutrino masses, 
and could provide the origin of the observed texture of the light neutrino mass matrix. Finally, we 
derive general formulae for sneutrino-antisneutrino oscillations and sneutrino flavor-oscillations. 
Unfortunately, neither oscillation phenomena is likely to be observable at future colliders. 



1 Introduction 



The Standard Model of particle physics provides a remarkable description of the funda- 
mental interactions of elementary particles at energy scales of order 100 GeV and below. 
Precision tests at LEP, the Tevatron and other lower energy colliders have detected no sig- 
nificant deviations from the predictions of observed electroweak phenomena [1]. Although 
the scalar sector responsible for electroweak symmetry breaking has not yet been discov- 
ered, the precision electroweak data is consistent with the Standard Model including a 
scalar Higgs boson of mass 114 GeV < rrih < 182 GeV at 95% CL. Despite its successes, 
the Standard Model is widely acknowledged to be only a low-energy effective theory, to 
be superseded (most likely at the TeV energy scale) by a more fundamental theory that 
can explain the puzzling large hierarchy between the energy scale that governs electroweak 
symmetry-breaking and the Planck scale [2]. 

Numerous proposals for a more fundamental theory that supersedes the Standard Model 
have been advanced over the last thirty years [3]. Low-energy supersymmetric theories (in 
which supersymmetry breaking effects of order the TeV scale are ultimately responsible for 
electroweak sjTumetry breaking) are perhaps the most well-studied framework for TeV-scale 
physics beyond the Standard Model [4-6]. The simplest supersymmetric extension consists 
of the particle content of the two-Higgs-doublet extension of the Standard Model and its 
supersymmetric partners. In addition to the supersymmetric interactions of the particle su- 
permultiplets, one adds the most general set of soft-supersymmetry-breaking terms, which 
parameterizes the unknown dynamics responsible for supersymmetry breaking [7,8]. The 
resulting minimal supersymmetric Standard Model (MSSM) yields a rich phenomenology 
of new superpartners and interactions, which if present in nature is poised for discovery at 
the Tevatron and/or Large Hadron Collider (LHC). 

Although no significant deviations from Standard Model predictions have been observed 
at colliders, there is of course one definitive set of observations that are in conflict with 
(the minimal version of) the Standard Model — the observation of neutrino mixing and its 
implications for neutrino masses [9]. Since neutrinos are strictly massless in the Standard 
Model, the latter must be modifled in order to incorporate the observed phenomena of 
neutrino oscillations. The simplest approach is to introduce a gauge invariant dimension- 
flve operator [lojll 

^5 = -^-f{ei,LlH,){euL'^H,) + H.c. , (1.1) 
where Hj is the complex Higgs doublet and Ll = {yi , £i) is the SU(2)-doublet of two- 
^Following refs. [7] and [6], we employ a convention where eyi = —1 = —£21- 
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component lepton fieldsjj where / and K label the three generations. 

After electroweak symmetry breaking, the neutral component of the doublet Higgs 
field acquires a vacuum expectation value, and a Majorana mass matrix for the neutri- 
nos is generated. The dimension-five term [eq. (11.11) ] is generated by new physics beyond 
the Standard Model at the scale A. Current bounds on light neutrino masses suggest 
that v'^/A < 1 eV [11, 12], or A > 10^^ GeV. A possible realization of eq. ffLTj) is based 
on the seesaw mechanism, which was independently discovered by a number of different 
authors [13,14]. In the seesaw extension of the Standard Model [14], one simply adds 
SU(2)xU(l) gauge singlet neutrino fields v'j^ and writes down the most general renormal- 
izable couplings of u'j^ to the Standard Model fields: 

^seesaw = -ei.vl' H.L^vj;' - \M''vl'vi' + H.c. (1.2) 

If ||M|| ^ f , then at energy scales below M a dimension-five operator of the form given by 
eq. (II. ip is generated. 

The MSSM is a minimal extension of the Standard Model. Nevertheless, there is a po- 
tential source for lepton-number violation and hence neutrino masses. Unlike the Standard 
Model, it is possible to construct renormalizable operators that violate lepton number and 
baryon number [15]. In their most generic forms, such operators would lead to extremely 
fast proton decay in confiict with the observations. The traditional solution is to intro- 
duce a discrete symmetry called R parity [16] that distinguishes Standard Model particles 
and their superpartners. In the R-parity-conserving (RPC) MSSM, neutrinos are massless 
just as in the Standard Model. Thus, one way to incorporate massive neutrinos in the 
RPC-MSSM is to formulate a minimal supersymmetric extension of the seesaw-extended 
Standard Model [17-21]. An alternative approach is to choose a different discrete symmetry 
that preserves baryon number but violates lepton number [22]. In such an R-parity- violating 
(RPV) MSSM, a Z3 baryon triality guarantees that baryon number is conserved by the 
renormalizable operators of the model (hence preventing fast proton decay). This approach 
has the advantage that no new fields beyond those of the MSSM need to be introduced. 
However, certain RPV (lepton-number- violating) couplings must be taken to be quite small 
in order to explain the scale of neutrino masses [23-25]. 

In this paper, we shall consider the minimal supersymmetric extension of the seesaw- 
extended Standard Model [17-21]. In this model, neutrino masses and mixing are governed 
by the same seesaw mechanism originally introduced into the (non-supersymmetric) Stan- 
dard Model. In the supersymmetry-extended model, new lepton-violating phenomena enter 

^To translate the two-component spinor product if into four-component spinor notation, 
see [Appendix A 
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due to additional effective lepton-violating operators generated by soft-supersymmetry- 
breaking. Such effects govern the behavior of the neutrino superpartners — the sneutrinos. 
Thus, the supersymmetric seesaw model provides new sources for lepton-number-violating 
phenomena. For example, sneutrinos and antisneutrinos can mix due to effective AL = 2 
operators [18,26]. Although such mixing effects are expected to be quite small, there are 
some scenarios in which sneutrino mixing phenomena could be observed in future collider 
experiments [18,27]. Sneutrino mixing also contributes a significant one-loop correction 
to neutrino masses and could be partially responsible for the observed pattern of neutrino 
masses and mixing [18,25,28]. The supersymmetric seesaw can also introduce lepton-flavor- 
violation and CP-violating effects due to the non-trivial flavor structure of the seesaw inter- 
actions [19,20,29]. Such phenomena are exhibited in the flavor oscillations of the charged 
sleptons [30] and the sneutrinos, respectively. Moreover, new one-loop processes contribute 
to £^ — £"^7 and electric dipole moments, and provide interesting constraints on the model 
parameters. 

In Section O we introduce the Lagrangian for the three-generation supersymmetric see- 
saw model, focusing on the interaction of the lepton and Higgs superfields. Our notation 



for fermion fields are described in Appendix A In Section [3l we derive the mass matrices 
for neutrinos and squared-mass matrices for the sneutrinos. In the limit of M ^ f , one can 
use perturbation theory to obtain accurate analytical expressions for the diagonalization of 
the effective mass and squared-mass matrices for the light and heavy neutral fermion and 
scalar states, respectively. The origin of a non-decoupling contribution to sneutrino masses 



noted in Section [3] is provided in Appendix B In Section |U we examine the constraints on 
the lepton-number conserving parameters of the model due to the observed g — 2 of the 
muon, the (unobserved) electric dipole moment of the electron, and the unobserved radia- 
tive decays of charged leptons. In Section El constraints on the lepton-number violating 
parameters of the model are obtained based on observed neutrino mass and mixing data. 
The general theory and phenomenology of sneutrino oscillations and mixing are addressed 
in Section [6l Our conclusions are given in Section [71 Although the neutrino are most easily 
treated as two-component spinor fields, it is convenient to present the Feynman rules of 



the model using four- component spinor notation. In Appendix A , we demonstrate how 
to translate between two-component and four- component spinor notation in the interac- 
tion Lagrangian. The relevant Feynman rules needed for the computations of this paper 
are listed in Appendix C Finally, some order of magnitude estimates for the contribu- 
tions to one-loop neutrino masses (relevant for the discussion of Section 15.11) are provided 
in [Appendix D 
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2 Lagrangian and the scalar potential 



In this section, we examine the terms of the Lagrangian that contribute to the masses and 
the non-gauge interactions of the neutrinos and sneutrinos. That is, we focus on terms 
that involve the charged leptons, neutrinos, charged sleptons, sneutrinos and the Higgs 
fields. The relevant superfields (denoted with hats above the corresponding field symbol) 
are specified in Table [H 



Table 1: 
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The electric charge (in units of e) is given by Q = T3 + y/2, where Y is the hypercharge 
specified above. The index j labels components of the SU(2) doublets with T3 = ±1/2 
for j = 1,2 respectively (and T3 = for the SU(2) singlets). The fermionic partners 
can be viewed either as two-component fermion fields or the left-handed projections of 



four-component fermion fields, as explained in Appendix A The index / = 1,2,3 labels 
three possible generations of charged lepton and neutrino superfields. The notation for the 
scalar field components of the hypercharge- zero superfield is motivated by the fact that in 
the lepton- number-conserving limit, R and possess the same lepton number (which is 
opposite in sign to that of L) . Consequently, ul and ur possess identical lepton numbers 
[cf. eq. (E3D]. 

The most general (renormalizable) form of the superpotential involving the lepton and 
Higgs superfields in the R-parity-conserving extended MSSM is given by: 

W = eij {fiHlH] - Y/-^HlL^jR^ + YI^H^L]N^) + iM^-^N^N^ , (2.1) 

where Yi and Y^ are complex 3x3 matrices, M is a complex symmetric 3x3 matrix and 
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/i is a complex parameterlfl In addition, there are soft-supersymmetry-breaking terms that 
involve the scalar field components of the above superfields. Before writing these terms 
explicitly, it is convenient to perform field redefinitions of the (charged and neutral) lepton 
superfields: 

V ^V[-^L\ ^V^-^R\ -^Vl/N\ (2.2) 

where V^, Vr and Vat are 3x3 unitary matrices. Note that the kinetic energy terms (and 
the couplings of the lepton superfields to the gauge fields) are invariant under the above 
unitary transformations. However, the coefficients of the terms of the superpotential are 
modified: 

Y.^VIY^Vr, Y^^VIY.Vn, M^V^MVm. (2.3) 
We shall choose Vl, Vr and Vat such that: 

Vf F,Vr = diag(Fe , Y^ , Yr) , (2.4) 
V^MVn = diag(Mi , M2 , Ms) , (2.5) 

where the elements of the two diagonal matrices above are real and non-negative. It is 
always possible to find unitary matrices Vl and Vr such that eq. (12.41) is satisfied — this 
is the singular value decomposition of an arbitrary complex matrix [31]. Likewise, it is 
always possible to find a unitary matrix Vn such that eq. (12.51) holds — this is the Takagi- 
diagonalization of an arbitrary complex symmetric matrix [31-33]. Thus, the redefinition 
of the lepton superfields [eq. (12. 2p ] implies that one can assume from the beginning without 
loss of generality that Yi and M are real non-negative diagonal matricesjj Note that the 
(transformed) Yy is in general an arbitrary complex 3x3 matrix. 

We next introduce the most general set of R-parity-conserving soft-supersymmetry 
(SUSY)-breaking terms (following the usual rules of [34]) involving the slepton, sneutrino 
and Higgs fields: 

VsoYT = mlHl*H] + m^^HrH^ + (mi)^-^Zf Zf + {miy-^R'*R' + (m^)^-^F*iV^ 

- {miy^N'N-^ + eij (ml^HlH] + A'/hIl'^R^ + A'J H^V^N^^ + H.c] , (2.6) 

where m|, m|, and m\ are hermitian matrices, m\ is a complex symmetric matrix and 
Ai and Ay are complex matrices. In general, these 3x3 matrices do not take a simplified 



With the convention for as specified in footnote [U it is convenient to insert an extra minus sign in 
front of Yi in eq. (|2.ip . This ensures that in a basis where is a real positive diagonal matrix, the charged 
lepton masses are also positive. Note that this convention differs from the one adopted in ref. [7]. 

^After electroweak symmetry breaking, eq. (j2.4p corresponds to working in a basis in which the charged 

lepton mass matrices are (real) non- negative and diagonal. 
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form in the basis defined by eqs. (12 ■4p and (I2.5p . The total scalar potential is made up 
of three contributions: the F-terms, which are derived from eq. (12. ip . the D-terms, which 
arise from the gauge interactions, and and the soft SUSY-breaking terms, which have been 
specified in eq. (12.60 . The total scalar potential is then given by: 



V = Vf + Vd + VsoFT , where Vp = ^ 



dW ^ 



(2.7) 



and the sum over i is taken over all scalar components of the corresponding superfields. 

The Yukawa couplings of the leptons and the Higgs fields and the corresponding fermion 
mass terms are derived from eq. (12.10 using the well-known formula [6,7]: 

'd^W\ 



_ (/> _ c/> _ 1 \p 

ij 



-ipiipj + H.c. 



(2.8) 



where the ipi are the two-component fermion field superpartners of the corresponding 0,, 
and W[(f)] is the superpotential function with superfields replaced by their scalar com- 
ponents. After electroweak symmetry breaking, the neutral Higgs fields acquire vacuum 
expectation values]^ 

(ff;> = ^, (*^'> = ^- (2.9) 

where v"^ = vf + V2 = (246 GeV)^ and tan/5 = V2/V1. Inserting the Higgs field vacuum 
expectation values into eqs. (12. 7p and (12.81) . one can isolate the terms of the Lagrangian 
that are quadratic in the scalar fields and fermion fields, respectively. These terms yield 
squared-mass matrices for the charged sleptons and sneutrinos and mass matrices for the 
charged leptons and neutrinos. In the basis defined by eq. (12.40 . the charged lepton mass 
matrix is diagonal, with diagonal elements m^i = viY/ / \/2. 

In general, the diagonalization of these mass matrices cannot be performed analytically, 
and one must resort to numerical techniques. However, the large hierarchy between neutrino 
masses and charged lepton masses strongly suggests that the parameters Mj ^ f , in which 
case an analytic perturbative diagonalization permits one to isolate the light (s)neutrino 
sector and integrate out the superheavy (s) neutrino sector, whose particle masses are of 
order the Mj. This procedure was carried out for the CP-conserving one-generation model 
in ref. [18]. In Section [3l we shall generalize this analysis to the most general (potentially 
CP-violating) three-generation model. 

First, we clarify the expected magnitudes of the parameters of the model: 



'''We define the overall phases of the neutral Higgs fields, Hi and such that the corresponding 
vacuum expectation values wi,2/V2 are real and positive. 
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1. We assume that the Yukawa couphngs Yj'^ satisfy 1^ 

linil <0(1). (2.10) 



The Majorana mass M is much heavier than the electroweak scale (seesaw mecha- 
nism [13]) 

||M||>w. (2.11) 



3. Ahhough yU is a supersymmetric parameter, we require it to be of a similar order to 
the low-energy supersymmetry-breaking scale, Msusy 



/i~MsusY. (2.12) 



4. The non-singlet soft SUSY-breaking squared-masses are of a similar order to the 
supersymmetry-breaking scale: 

\\<\\ ~ ~ MlusY • (2.13) 



5. The parameters and A^, are unconnected to electroweak symmetry breaking at 
tree-level. However, these parameters generate a mass-splitting between sneutrinos 
and antisneutrinos. The latter contributes via loop corrections to neutrino mass 
splittings, which are experimentally constrained. One expects that [36]: 

P.ll < Msusy , < Msusy||M|| , (2.14) 

although these parameters could conceivably be larger by as much as a factor of 
10^ [18]. Large A^, also leads also to large corrections to charged slepton masses. 
Thus, to avoid unnatural fine-tuning in order to prevent charged slepton masses from 
being larger than about 1 TeV, one again expects that A^ cannot be much larger 
than the supersymmetry-breaking scale. The impact of the one-loop effects of on 
charged lepton radiative decays and the Higgs mass parameters also yield constraints 
and imply that the bound on given by eq. fl2.14p cannot be significantly relaxed. 

6. The singlet soft SUSY-breaking parameter is also unconnected to electroweak 
symmetry breaking at tree-level. However, the one-loop corrections to the Higgs 
mass parameters depend quadratically on m^, so to avoid unnatural fine-tuning of 

^The Euclidean matrix norm is defined by \\A\\ = [tr(A'l'A)] = 
matrix elements are given by aij . 



E 



1/2 



for a matrix A whose 
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the electroweak symmetry breaking scale, one expects that cannot be much larger 



\2 



than (1 TeV) . This expectation is confirmed in Appendix B in which case 



IKII ^ Ms^usY • (2.15) 

If significant fine-tuning of the electroweak scale is allowed (as in the split-super- 
symmetry [37] approach), then the constraints on are significantly relaxed. The 
one-loop effects of mj^ on physical observables are rather mild, even as ap- 
proaches ||M^||. For example, in ref. [38], the one- loop corrections to Higgs masses 
in the seesaw-extended MSSM are found to be large and negative if \\mj^\\ , ||7Ti^|| ~ 
||M^||. However, these corrections become negligible once these soft-SUSY-breaking 
masses are taken somewhat below the seesaw scale. 

Thus, we shall present results in this paper that allow for the possibility that: 

llm^ll ~ IIM^II . (2.16) 

If eq. (I2.16P holds, then remnants of the heavy neutrino/sneutrino sector can survive 
in the effective theory of the light sneutrinos. The origin of this non-decoupling effect 



is explored in Appendix B 



Although naturalness demands that the scale of low-energy supersymmetry-breaking, 
^susY, should be (roughly) of 0{v), the absence of observed supersymmetric phenomena 
(and a light CP-even Higgs boson) suggest that MgusY may be somewhat larger, of order 
1 TeV. Nevertheless, in eqs. fl2.12p - fl2.15p . one could substitute Msusy with v; the results 
of this paper are consistent with either choice. 



3 The (s) neutrino (squared-) mass matrices 

In this section, we examine in detail the neutrino mass matrix and the sneutrino squared- 
mass matrix. In a three-generation model, the neutrino mass matrix is a 6 x 6 complex 
symmetric matrix, which can be written in block (partitioned) form in terms of 3 x 3 
matrix blocks. The sneutrino squared-mass matrix is a 12 x 12 hermitian matrix, which 
can be written in block (partitioned) form in terms of 6 x 6 matrix blocks. Each of these 
6x6 matrices can be further partitioned in terms of 3 x 3 matrix blocks. In order to 
accommodate the proliferation of matrices of dimension 3, 6 and 12, we adopt a notational 
device that allows the reader to instantly discern the dimension of a given matrix. Thus, 
we use a boldface capital letter {M) to denote a 12 x 12 matrix, a calligraphic letter (A^) 
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to denote a 6 x 6 matrix, and a Latin letter (M or m) to denote a 3 x 3 matrix. Latin 
letters will also be used to denote (scalar) mass parameters, with appropriate identifying 
subscript or superscript labels to distinguish these from the 3x3 matrices introduced in 
Sections [2] and [31 Following the conventions of Section [2], we shall employ subscript and 
superscript upper case Latin indices /, J, K as generation labels that run from 1 to 3. 
Lower case Latin indices i, j, k are employed for other purposes, either as SU(2) gauge 
indices or as labels representing the six light sneutrino mass eigenstates. Other subscripts 
appearing in this section will be used to distinguish among different matrix quantities. 



3.1 The neutrino mass matrices 

Working in a basis where M is a diagonal matrix [cf. eq. (12. Sp ], we begin by analyzing the 
neutrino mass matrix. The resulting terms quadratic in the neutrino fields are given in 
terms of two-component fermion fields l| by: 

- = i {v2V2Yj'ulul' + M^'uful' + H.c.) = i {ul uf) (^"''^^ + H.c. 

(3.1) 

The neutrino mass matrix A^j^ is a 6 x 6 complex symmetric matrix given in block form 
by: 

[ mo A 
M,^\ , (3.2) 

y M J 

where the 3x3 complex matrix 

mn = V2Y,/V2 (3.3) 

generalizes the neutrino Dirac mass term of the one-generation model [cf. eq. ( lA.SIl . 

Provided that ||M|| ^ HfnoW [as suggested by eq. (12.111) ]. J^i, is of a seesaw type [13]. 
The neutrino mass matrix can be Takagi block-diagonalized [21,25,33] as follows. Introduce 
the 6x6 (approximate) unitary matrix: 

U=\ ' I , (3.4) 

where 1 is the 3x3 identity matrix. 



^ In [Appendix Aj we show how to rewrite eq. (j3.ip in terms of four-component neutrino fields. However, 
the two-component formalism is more economical, so we adopt this notation in what follows. 
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One can check that: 



, 1 + 0{m%M-^) 

1 + C(m|)M-^) 



We define transformed (hght and heavy) neutrino states vi and z/^ by: 




(3.6) 



By straightforward matrix multiphcation, one can verify that 

\ 0{m%M~^) M+l{M-^m^^mD + m'^m*jjM-^) + 0{m%M- 

(3.7) 

At this stage, we can identify an effective (complex symmetric) mass matrix M^,^ for 
the three hght (left-handed) neutrinos with respect to the {i/^j-basis: 

M^^ ~ -moM-^ml . (3.8) 

To identify the physical light neutrino states, we must perform a Takagi-diagonalization 
of My^. This is accomplished by introducing the unitary MNS matrix [39], Umns, via 



where the (z//)P^y'^ [J = 1, 2, 3] denote the physical light neutrino fields. Umns is determined 
by the Takagi-diagonalization of M,y^: 

f^MNS^i^^f^MNS = diag(m^^^ , m^^^ , m^^^) , (3.10) 

where the 777/ J cLFG the (real non-negative) masses of the light neutrino mass eigenstates. 

For completeness, we examine the effective mass matrix of the heavy neutrino states. 
Although M is diagonal by assumption, the lower right-handed block in eq. (13. 7p is no 
longer diagonal due to the second-order perturbative correction. However, we do not have to 
perform another Takagi-diagonalization, since the off-diagonal elements are of 0{mjyM~^), 
and would only affect the physical (diagonal) masses at order 0{mjyM~^), which we neglect. 
The corresponding mixing angles would be of (9(m|,M~^), which we also neglect here. 
Thus, we identify the physical heavy neutrino mass eigenstates to leading order by: 

(^./)phys ^ ^c/ ^ (3.11) 
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with masses 

'fl + li^EK'l') . (3-12) 



where the Mi are the diagonal elements of M in our chosen basis. 



'I J 



3.2 The sneutrino squared- mass matrices 

We now turn to the sneutrino sector. It is convenient to separate out various pieces that 
comprise the F-term contributions to the scalar potential [eq. ( 12. 7p ]: 

Vp = K + ^M + Kther, (3.13) 

where K = Z]i=L{ \dW/d(j)i\^ and = \dW/dHl\'^ ultimately contribute to the sneu- 
trino squared-mass matrix, whereas Kther (which involves derivatives of the superpotential 
with respect to the other scalar fields) makes no contributions to tree-level sneutrino masses. 

As a pedagogical exercise, we first analyze the supersymmetric limit. Although super- 
symmetry-breaking is required in the MSSM to generate electroweak symmetry breaking, 
one often finds supersymmetric-like relations between the fermion and sfermion sectors in 
the limit of vi = V2 and /i = 0, i.e. for = Vd = 0. Thus, in the following computation 
the supersymmetric limit corresponds to taking the total scalar potential [eq. (12. 7p ] to 
he V = Vy. To analyze the contributions of Vy to sneutrino masses, we can employ the 
following trick. Focus on the following two terms of the superpotential: 

Consistent with eq. (13.61) . we redefine the neutrino superfields as follows: 

(3.15) 







(Lu\ 















where the unitary matrix U is given by eq. (13. 4p . Defining the matrix H = i/|yj/, the effect 
of eq. (I3.15P is to transform intccl 

Wy ~ l{HM-^H^y-^L{^L(^+l [M" + \{M-^H^H + H^H*M-y^] NInI+0{H^M-^) , 

(3.16) 



^Strictly speaking, this is not a permissible transformation, since W must be holomorphic in the super- 
fields, whereas eq. p.l6p is a function of both i7| and However, since we ultimately set 77| — V2/V^ 
and only take derivatives of Wv with respect to Lu and Nh, the procedure outlined here yields correct 
results. 
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where there is an imphcit sum over / and J. In deriving eq. (13.161) . we have used the fact 
that M^'-' is a non- negative diagonal matrix. Setting if| = V2/\/2 and using eq. (12. 7p . we 
can directly make use of eq. (13.161) to isolate the contributions to the sneutrino squared-mass 
matrix that arise from Vy. 

- ^^ass = L\,MI,Lu + NlMl^Nh , (3.17) 

where the 3x3 hermitian matrices M'j-^^ and M^|^ are given by: 

Ml^ = m*ijM-^m^^mDM-^ml + 0{m%M-^) , (3.18) 
M^t^ = m^^mo + |(Mm^m^M-^ + M-^m]^m*jjM) + 0(m^M-2) . (3.19) 

Moreover, the effective light and heavy neutrino mass matrices, My^ and My^, can also be 
derived by inserting eq. (13.161) into eq. (12. 8p . As expected, the resulting neutrino mass ma- 
trices are related in a supersymmetric way to the sneutrino squared-mass matrices obtained 
in eqs. flXTBD and (ITOD : 

Ml, = MlMy,, Ml,^ = MlMy^. (3.20) 
In particular, in the supersymmetric limit, 

t^MNS f^MNS = diag {ml^^ , ml^^ , mlj , (3.21) 

which implies that the light neutrino and sneutrino masses coincide. 

We now turn to the complete calculation of the sneutrino mass matrix. Although one 
could perform the computation with respect to the basis of sneutrino states defined by 
eq. (I3.15p . this basis is not especially convenient. This is due to the fact that the effec- 
tive squared-mass matrix of the light sneutrinos is dominated by supersymmetry-breaking 
effects. In particular, the supersymmetric contribution of 0{m\)M~'^) [cf. eq. (I3.18p ] is 
completely negligible relative to the supersymmetry-breaking contributions. Thus, there is 
no advantage to performing in the sneutrino sector the same change of basis used to isolate 
the effective mass matrix of the light neutrinos. Hence we will write the 12 x 12 hermitian 
sneutrino squared-mass matrix in block form as: 

= i( 4 II I . (3.22) 



2 

NN 




where 0l = {Li , L\)^ and (pN = {N , A^*)^ are six-dimensional vectors. The 6x6 hermitian 
matrices M.\li M.%n ^^"^ 6x6 complex matrix A^|^ can be written in block partitioned 
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= m\ + \Ml cos 2(3 + m})m^ 




= + m% + m}j^mD , 


Km 


= m^M, 




= -X^mn, 


^NTN 


= -2ml , 




= 0, 



form as: 

Ml, . ( I , (3.23) 

where the subscripts A and B can take on possible values L and N [this labeling allows 
one to keep track of the origin of the various matrix blocks]. The M^^^ are 3x3 hermitian 
matrices and the M'^t^ 3x3 complex symmetric matrices, for A = L , N . There are 
no restrictions on the 3x3 complex matrices Af^t_B ^"^^ ^atb A ^ B. 

Adding up the contributions of V^, Vd and Vsoft to the sneutrino masses yields: 

(3.24) 
(3.25) 
(3.26) 
(3.27) 
(3.28) 
(3.29) 

where we have introduced the complex 3x3 matrix parameter Xj, by the following definition: 

X.niD = ^ {V2A, + fi*vM . (3.30) 
v2 

A quick check of the supersymmetric limit confirms the expected relation between the 
neutrino mass matrix and the sneutrino squared-mass matrix: 

, m*nnv^ m*nM 
MlM. = I ^ ^ ^ I . (3.31) 

As noted above, because of the dominance of supersymmetry-breaking contributions to the 
light sneutrino masses, the diagonalization of the light neutrino mass matrix and the light 
sneutrino squared-mass matrix are completely independent. 

Under the assumptions of eqs. fl2.10l) -( l2.15l) . the 12 x 12 sneutrino mass matrix, written 
in terms of 6 x 6 matrix blocks with estimated magnitudes, 

MNI ^" (3.32) 
(A<i„)* M%„J \0(vM) 0(M^' 

also exhibits a seesaw type behavior, analogous to the seesaw type mass matrix [eq. (13.21) 
of the neutrino sector. Following the standard procedure for diagonalizing such matrices 
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see ref. [25]), we introduce a 12 x 12 unitary matrix: 

(3.33) 



2-'^'-LN'^'"-NNV'-LN ) ■^'"■LN-""-NN 

where X is the 6x6 identity matrix. One can easily compute: 

Hence, the effective 6x6 hermitian squared-mass matrix for the hght sneutrinos reads: 

Ml ^ Ml^ - MIj,M-A {Mlj^y + 0{v'M-') , (3.35) 

analogous to the light effective neutrino mass matrix of eq. (13.81) . Likewise, the effective 
6x6 hermitian squared-mass matrix for the superheavy sneutrinos reads: 

M\ ^ Ml^ + \ [M-^\[Ml^)^Ml^ + {Ml^)^Ml^M-^\] + 0{v'M-') , (3.36) 

where for completeness, we have exhibited the 0{v'^) corrections to the leading term. As 
expected, the masses of half of the sneutrino eigenstates are of order the electroweak sym- 
metry breaking scale, whereas the other half are superheavy, of order M. 

Following the notation of Table [1], the (complex) sneutrino interaction eigenstates are 
denoted by: vl = Li and = A^*. The latter convention reflects the fact that in the 
lepton-number conserving limit of M^'^ = m\ = 0, the lepton numbers of Ul and are 
identical, as previously noted. (Of course, the limit of interest in this paper, ||M|| ^ v, 
is very far from the lepton-number conserving limit.) In analogy to Ui and Uh, we define 
transformed (light and heavy) sneutrino states Ui and Uh by: 

(3.37) 

where = (z/^ , z/|)'^ and (f)h = (j^l , i^h)'^ are six-dimensional vectors. Sneutrino- 
antisneutrino oscillations are a consequence of the AL = 2 elements in the light and heavy 
sneutrino squared-mass matrices -M^ and -M^, and are governed by M^^j^ and M^^^ 
(note that M^^^, which would also violate lepton number by two units, is zero). 

Using the form of M.\b B = L ot N) given by eq. (I3.23P with the given in 
eqs. (I3.24p -( !3.29p . the effective 6x6 hermitian squared-mass matrix for the light sneutrinos 
[eq. fl3.35p ] is given by: 

Ml ^1 ^ I , (3.38) 

Miy {Mi^y 
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where the lepton-number-conserving (LC) and lepton-number-violating (LV) matrix ele- 
ments are given by: 



Mlc = ml + \Ml cos 2(3 + m^m^ - m^M(M2 + m^Y^Mm^^ + 0{v'^M~'^) , (3.39) 

Mly = moMiM^ + m^*)-^m5xj + X^nioiM^ + m%)-^MmJ) 

-2mDM{M^ + m\*)-^m\{M'^ + mly^Mml + O^v^M''^) , (3.40) 

under the assumption that and can be as large as indicated in eqs. fl2.14p and 
fl2.16p . Note that M|(^ is a 3 x 3 hermitian matrix, and is a 3 x 3 complex symmetric 
matrix. Moreover, although M is a diagonal matrix with real positive entries [cf. eq. (12.51) ]. 
m\ can be any 3x3 hermitian matrix, not necessarily diagonal nor real. The M ^ oo 
limit of eqs. fl3.39p and fl3.40p is noteworthy. In this limit, = and the lepton- 

number-violating effects completely decouple, as expected. If in addition = 0, then 
M|p = m\ + |M| cos 2/5, which reproduces the well known 3x3 light sneutrino squared- 
mass matrix of the MSSM. However, according to eq. fl2.15p . m^M~^ ~ C^(l) is possible, 
in which case deviates from its MSSM value by a quantity of 0{v'^) even in the exact 
decoupling limit of M —>■ oo. The origin of this non-decoupling behavior is explained in 



Appendix B As a result of this non-decoupling phenomenon, remnants of the heavy sector 
of the seesaw mechanism may survive in the effective theory of light sneutrinos. These 
non-decoupling effects can be detected in principle through measurements of the sneutrino 
and charged slepton properties. 

The physical light sneutrino states can be identified by diagonalizing Ai^^. Note that 
if M|y = 0, then the eigenvalueqj of Ml^, are doubly degenerate, corresponding to the fact 
that the conserved lepton number implies that the six light sneutrino states are comprised 
of three sneutrino antisneutrino pairs. If M|y 7^ 0, then lepton number is violated and the 
sneutrinos and antisneutrinos can mix. This mixing splits the degenerate pairs and yields 
(in general) six non-degenerate light sneutrinos. In particular, the resulting sneutrino 
mass-eigenstates are self-conjugate real fields, which we denote by , 5*2 , . . . , 5*6. 

To determine the 5*^ in terms of the interaction sneutrino eigenstates, one must compute 
the 6x6 unitary matrix W that diagonalizes A4?^: 

W^MlW = diag {ml ,ml,..., mjj . (3.41) 

Noting that STW^S = Ml*, where S = (§ J), it follows that if W satisfies eq. (ICTj) 
then so does SW*. However, the unitary matrix that diagonalizes is unique up to 



^Under the assumption that R-parity is not spontaneously broken, the (real) eigenvalues of the hermitian 
matrix are non-negative. 
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a multiplication on the right by a unitary matrix Ud that is arbitrary within a subspace 
of degenerate eigenvalues and is otherwise diagonal. Denote the set of all such unitary 
matrices by S. Hence, one can conclude that EW* = WUd for some Ud € S. Since W is 
unitary, Uo = W^SW*, and it follows that UdW^ = I. That is, Ud must be a symmetric 







1 /2 

unitary matrix. It then follows that the matrix W = Wf/^ satisfies W = SW" 

Thus, without loss of generality, we may drop the primed superscripts and impose the 
constraint W = SW* on the diagonalizing matrix that satisfies eq. (13.411) . It then follows 
that W has the following form: 

fx \ 

W = , (3.42) 

y X* -lY* ) 

where X and Y are 3x3 complex matrices that satisfy: 

XX^ + YY^ = 1 , XX^ = YY^ , (3.43) 

Re{X^X) = Re{Y^Y) = | , lm{X^Y) = , (3.44) 

due to the unitarity of W. Consequently, the relation between the sneutrino interaction- 
eigenstate fields z// and the six self-conjugate sneutrino mass-eigenstate fields is given 
by: 

6 3 

^'i=Yl ^''^'^ = Yl {^'""Sk + tY'"" Sk+s) , (1 = 1, 2, 3) . (3.45) 

k=l K=l 

One can then invert eq. (13.451) [using eqs. (I3.43P and (13.441) ] to obtain: 

3 3 
7=1 7=1 

(3.46) 

Indeed, the Sk are self-conjugate real fields as noted above. 

Since ~ O^v"^) and M|y ~ 0{v^M~^), the mass-splittings of the would-be 

sneutrino-antisneutrino pairs are expected to be very small, of order a typical neutrino 
mass. To compute the magnitude of the corresponding mass-splittings, we can employ per- 
turbative techniques to evaluate the eigenvalues of A^?^ [eq. (13.381) ]. First, we diagonahze 
the sub-matrix M|^: 

QlMlcQo = D = dmg{d^ , , d^) , (3.47) 



""^^We define U^'^ G iS to be the unique square root of Ud that is symmetric and unitary. This is 
accomplished by noting that there exists a (unique) real symmetric matrix Ti. such that Ud ~ exp(i7i). 

1 /2 

Then, Uj^ = cxp{iH/2). Note that there is still some freedom left in the choice of W' , which is unique up 
to a multiplication on the right by a real orthogonal matrix that is arbitrary within a degenerate subspace 
and is otherwise diagonal. 
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where Qo is a 3 x 3 unitary matrix, and the eigenvalues dj are reaL Note that Qo is not 
unique. In Section 14.31 we will argue that the bounds on the radiative flavor- changing 
charged lepton decay i"^ — > i^'j imply that matrix is very close to a diagonal form. In 
the limit of diagonal M|^, we shall take Qq = 1. We can then determine the off-diagonal 
elements of Qq by writing ~ diag(mj , , rrig) + where m'j^^ is a matrix made 
up of the off-diagonal elements of M^q, and Qo — 1 + Q'o, where = —go- By assumption, 
the matrix elements of m'j^Q are much smaller than the mj, and the matrix elements of go 
are much smaller than unity. Thus treating eq. (13.471) to flrst order in the small quantities, 
we can solve for the off-diagonal elements of go in terms of the elements of w^ic* ^"^^ 
m|. Since at flrst order mj = dj, it follows that: 

(Qo)/.-^^, ly^J. (3.48) 

The diagonal elements of Qo can then be determined to the same order by using the 
unitarity of Qq. In the remainder of this section, we will not make any assumption regarding 
the size of the off-diagonal elements of M|^, in which case eq. (13.481) does not apply and 
Qo must be obtained numerically from eq. (I3.47p . 

In the following, it will be convenient to deflne 

Q = QoT (3.49) 

where T is a 3 x 3 diagonal matrix of phases given by: 

T ^ diag , e-""-/' , e-^^^^/^) , 0j ^ arg {Q^MlyQo) . (3.50) 

Note that the right hand side of eq. (I3.47P is unchanged when Qq QqT, so that the 
unitary matrix Q can also be used to diagonalize M|^. It then follows that: 



D B*\ \ Mlc (MlyY 



B D I \ J \ Mly (M2 



LCJ 




(3.51) 



where B is the 3x3 complex symmetric matrix 

B = Q^MlyQ . (3.52) 

Due to the rephasing of Qq as specifled by eqs. (13.491) and (I3.50p . the diagonal elements 
of B are real and non-negative: Bjj = \Bjj\. This is the motivation for our choice of Q 
in the diagonalization of Ml^. Note that if is approximately diagonal, then Qo — 1, 
in which case (pj ~ arg[(M|y) jj]. Thus, unless the diagonal elements of M|y are non- 
negative, Q ~ T 7^ 1 in this limiting case. 
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Even though D ~ 0{v'^) and B ~ 0{v^M~^), the unitary matrix that diagonahzes & 
is not close to the identity matrix, due to the double degeneracy of the diagonal elements. 
In order to perform a perturbative diagonalization of we first introduce the following 
6x6 unitary matrix P, expressed in block form as: 



(3.53) 



A straightforward computation yields: 

, D + ReB -ImB , 

r^^V =1 I , (3.54) 

-ImB D-ReB 

which is a 6 X 6 real symmetric matrix. 

If the elements of the diagonal matrix D are non-degenerat such that dj — dj ~ 0{v'^) 
for all I ^ J, then the matrix V'^ &V can be diagonalized by a real orthogonal matrix TZ 
that is close to the identity: 

/ 1 + Rei? Imi? \ ^ , 

7^ = + 0{v^M-^) , (3.55) 

y ImR 1 - ReR ) 
where the 3x3 complex antisymmetric matrix R is of order (9(fM~^): 

Rij = -Rjj = , {I^J). (3.56) 
a J — aj 

One can check that: 

n^V^^VTZ = diag(m|^ , , . . . , ,m|J + 0{v*M-^) , (3.57) 

where the squared-masses of the light sneutrinos are given by: 

, 5.+3 = dj±\Bjj\ + 0{v^M~') , ( J = 1, 2, 3) , (3.58) 

and > rn'^j_^^. Note that the perturbations due to the off-diagonal elements of B 
contribute only to the (9(f^M~^) terms of the squared-masses. 



^^In general, we would expect the dj (which are the eigenvalues of M^q) to be non-degenerate. Even if 
the parameters and were proportional to the identity matrix at the high energy scale due to some 
flavor symmetry, this latter symmetry would not be respected by the corresponding low-energy parameters, 
due to flavor-violating effects that enter the renormalization group running. Moreover, the matrix mo is 
likely to reflect some of the flavor-violating effects of the model. Hence, any (near) degeneracy among the 
dj would be purely accidental. 
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Combining the results of eqs. fl3.5ip . fl3.53p and (13.551) . the hght sneutrino mixing matrix 
[defined in eq. (I3.4ip ] is given by: 




(3.59) 



■.Q{1- R) +0{v^M-^) . (3.60) 



Q{l + R) iQ{l-R) 
Q*{l + R*) -iQ*{l-R*] 

Comparing with eq. (13.420 . we identify: 

X = ^Q{l + R)+0{v^M-^), and Y = 

V 2 V 2 

Inserting these results into eqs. (13.450 and (I3.46P yields the desired (approximate) relations 
between the sneutrino mass eigenstates Sk and the interaction eigenstates u^. 

For completeness, we briefiy examine the modifications to eq. (I3.58P if some of the dj are 
degenerate. In this case, the diagonalizing matrix TZ is not close to the identity matrix, and 
the perturbative analysis above fails. Consider the case of dj = dj ^ dx, where {J, J, K} 
is some permutation of {1,2,3}. The first order shift in the eigenvalues of ^ will depend 
on Bjj as well as on the diagonal elements of B. However, the perturbations due to Bjk 
and BjK will only generate second-order shifts to the eigenvalues, which we neglect here. 
Thus, it is sufficient to solve the characteristic equation of ^ in the limit of dj = dj and 
Bjx = BjK = 0. In this limit, the characteristic polynomial factors into a product of two 
simpler polynomial factors!^ 



(A-rf 



K 



\B 



KK\ 



(A - di)' - (A - dj)' [\Bu\' + \Bjj\' + 2\Bjj\'] + \B 



2 

IJ 



BiiBjj\ 
(3.6 



The resulting sneutrino squared-masses are: 





±] 


.2 


\Bji 




\Bjj 


^ + 2\Bjj\'' + ^ 


/a" 


1 








\Bji 




\Bjj 


^ + 2\Bjj\^-^ 


/a" 


1 



1/2 

1/2 



^Sk,Sk+3 - dk±\BKK\ 



where 



A 



\B 



11 \ 



+ \B 



2\B 



iJ\ 



4 \ Bjj - BiiBjA 



(3.62) 

(3.63) 
(3.64) 

(3.65) 



The corresponding mixing matrix can be obtained by performing an exact diagonalization 
within the two-dimensional degenerate subspace, although we shall omit the details. 



-"^^In the case of a near degeneracy where dj ~ dj < ©(uAf"^), the quartic polynomial factor of the 
characteristic equation of S! contains a term linear in A — ^(rf/ + dj). In this case, the resulting expressions 
for m?g^ g^^^ and Sj^3 ^'"^ significantly more complicated than those presented in eqs. (|3.62p and 
(IM) . 
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Finally, in the very unlikely scenario where di = d2 = ds = d, all of the matrix elements 
of B contribute to the first order shifts of the eigenvalues of ^. To determine these shifts, 
put X = d + X in the characteristic equation of ^ to obtain a sixth order polynomial in x. 
No further perturbative simplification is possible, since all the terms of this polynomial are 
of the same order of magnitude. 

As expected, the mass-splittings of the would-be sneutrino-antisneutrino pairs are 
nonzero due to the presence of the lepton-number violating matrix M|y [cf eq. f l3.52p ]. If 



we denote the three sneutrino mass-splittings by (Amj^Jj 
then in the non-degenerate case. 



\B 



'di 



msj-msj^^ \ (for J = 1,2,3), 



(3.66) 



In the case of degenerate dj, the mass-splittings {Arrij^Jj also depend on the non-diagonal 
elements of B. 

It is instructive to examine the above results in a simplified one generation model. In 
this case, D = MIq and B = M^y are just numbers. In particular, m% is a real parameter 
and A^?^ is a 2 X 2 hermitian matrix, with eigenvalues 



m 



Sl,S2 



Mlc 



± \Mly\ 



ml + \Ml cos 2/? + 



M2 + 



± 



TV 



2|mBpM 
M2 + 



N 



M2 + 



m 



N 



(3.67) 



The corresponding sneutrino mass-splitting, Am^^ = \^S2 ~ ^sA, is given by 



Am 



2M^ 



m 



m^,{M^ + mlj) 



Mm% 



M2 + ml 



N 



(3.68) 



where m^^, = \mD\'^/M is the mass of the light neutrino and m^^ = ^{ms^ + ms^) is the 
average light sneutrino mass. If -C M, then eq. fl3.68p coincides with the result given 
in ref. [18] after taking into account a slight difference in notation]^ 

Assuming that ~ 0{vM), it follows that both terms on the right hand side of 
eq. (13.681) are of the same order, which implies that Am/>^ ~ 0{m^^). However, as noted 
below eq. (I2.14p . it is possible that could be as much as a factor of 10^ larger than its 
naive estimate [18], in which case the sneutrino-antisneutrino mass splitting could be three 
orders of magnitude larger than the corresponding light neutrino mass0 



-"^•^If we put = —M _Bjv and change the sign of (with the corresponding change in X^, [cf. eq. ()3.30|) ]). 
we recover the resuhs of ref. [18]. 

""^■^A similarly enhanced sneutrino-antisneutrino mass splitting also arises in the supersymmetric triplet 
seesaw model of ref. [40] . 
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The same set of manipulations described above can be carried out to obtain the cor- 
responding results for the effective 6x6 hermitian squared-mass matrix for the heavy 
sneutrinos [eq. (13.361) ]: 

Ml^l ^ ' + Oiv'M-') , (3.69) 

V -2ml mr j 

where the 3x3 hermitian matrix M'jj is defined by: 

Mjj = M^+m%+m^^mD + l{M^+m%)-^Mmlm*j^M+lMm^m}jM{M^+m%)-^ . (3.70) 

The physical heavy sneutrino mass-eigenstates are determined by diagonalizing A^|^. At 
leading order, the mass-eigenstates are mass-degenerate sneutrino/antisneutrino pairs, with 
masses and mixing angles (with respect to the basis in which M is diagonal) determined 
by the diagonalization of m^. The lepton-number violating off-block-diagonal matrix 
generates sneutrino-antisneutrino mixing, and yields mass-splittings between nearly degen- 
erate heavy sneutrino pairs of order Amj>^ ~ (9(m^M~^). 

The complex elements of the sneutrino squared-mass matrix govern CP-violating sneu- 
trino phenomena, due to the non-degeneracy of masses of the real and imaginary parts of 
the sneutrino fields. Following the discussion of the CP-properties of the sneutrino fields 
in Section [6], we find it convenient to define a new basis of sneutrino interaction eigenstates 
of definite CP. That is, we decompose the complex sneutrino fields into real and imaginary 
parts: 



''^7=2 



(3.71) 
(3.72) 



where the — ] superscripts indicate that the corresponding sneutrino eigenstates are 
CP-even and CP-odd. With respect to the CP-basis, 



(3.73) 

where V is the 6x6 unitary matrix introduced in eq. (13.531) . 

That is, with respect to the CP-basis, the effective squared-mass matrix for the light 
sneutrinos is given by: 

— 2 . . f ^<Mlc + Mly) -lm{Mlc + Mly) , 

Im(M|^ - Mly) Re{Mlc - Mly) 
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This is a real symmetric matrix (which is easily checked by recaUing that M^^ and 

are, respectively, hermitian and complex symmetric matrices), as the CP-basis consists of 

real self-conjugate scalar fields. 

If ImM|^ = ImM|y = 0, then the sneutrino mass-eigenstates are also definite eigen- 
states of CP. If in addition ReM|y 7^ 0, then the would-be sneutrino-antisneutrino pairs 
are organized into CP-even/CP-odd pairs of nearly degenerate sneutrinos [18]. 

Since is real symmetric, it can be diagonalized by a 6 x 6 real orthogonal matrix, 
Zc, via: 



(3.75) 



and the corresponding physical sneutrino mass eigenstates, 5*^ (fc = 1, . . . , 6), can be iden- 
tified as linear combinations of the CP-even and the CP-odd sneutrino eigenstates: 



(3.76) 



Matching with the notation employed by our discussion of sneutrino oscillations in Sec- 
tion [6l we note that the sneutrino interaction eigenstates, z/^, can be expressed in terms of 
the physical (self-conjugate) sneutrino mass eigenstates via: 



Comparing eqs. (13.451) and (13.771) . we can identify: 



which can be inverted to obtain: 



(3.77) 



[/,ir=l,2,3) 
(3.78) 



ZJ, = ^p2 



Re X -Im Y 



Im X 



Re F, 



(3.79) 



One can easily verify that the orthogonality of Zi, implies the unitarity of W defined in 
eq. (13.421) [and vice versa]. In particular, eqs. (I3.4ip and (13.751) imply that Zy = V^'W, in 
which case 

'^'^ JjW = >VtW=J, (3.80) 



ZTZy = W^V*V^W = W 
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after using the explicit forms for W and V. 

In summary, we have derived the hght effective sneutrino squared-mass matrix by ex- 
ploiting the seesaw mechanism in the sneutrino as well as in the neutrino sector. Our 
calculation is quite general under the parameter assumptions specified by eqs. (12.101) - 
(I2.15p . We found that Aif^ depends on two 3x3 matrix blocks, and M^y, given 

by eqs. (13.391) and (13.401) . respectively. In particular, M|y is responsible for the splitting 
of the masses of would-be sneutrino-antisneutrino pairs, or equivalently the mass-splitting 
of CP-even/CP-odd sneutrino pairs, t'^^^'*, in the CP-conserving limit. As we shall see in 
Sections H] and [5], the matrices M^^ and M^y provide a convenient parameterization for a 
number of interesting physical observables, such as neutrino masses and radiative lepton 
decays. 

4 Constraints on lepton number conserving parame- 
ters 

The input parameters that govern sneutrino mixing phenomena and sneutrino decays are 
encoded in matrices M|y and M^^^ given by eqs. (13.401) and (I3.39p . respectively [or, alterna- 
tively, in the physical sneutrino masses and the orthogonal matrix Zc, defined in eq. (13.751) . 
At present, apart from neutrino oscillations, only lepton number conserving processes are 
observed in current experiments. These processes constrain the entries of the lepton num- 
ber conserving matrix Ml^. In this Section we investigate bounds on the structure of 
imposed by the measurements of the muon magnetic moment anomaly, the — the elec- 
tric dipole moment (EDM) of the electron and the radiative flavor changing charged lepton 
decays, £^7. The latter have also been worked out in detail in ref. [21]. Additional 

constraints due to l~j ijiji'j' decays and /i-e conversion in nuclei are also relevant and 
have been analyzed in Ref. [21,41]. These constraints can yield further restrictions on the 
structure of MIq, although we shall not present this analysis here. 

We briefly summarize the constraints from current experiments relevant for the compu- 
tations presented in this Section. The most recent experimental measurement of the muon 
anomalous magnetic moment {a'^^^) exhibits a slight discrepancy [42] relative to the pre- 
dicted value of the Standard Model (ajf ). A recent theoretical review of the computation 
of the Standard Model prediction [43] yielded = a^^^ - = (2.94 ± 0.89) x 10-^ 
where all theoretical and experimental errors are added in quadrature, corresponding to 
a 3.3(7 effect. Thus, we roughly expect that the contribution to the muon anomalous 
magnetic moment from new physics beyond the Standard Model to be no larger than 
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iJ 1 f 2 £1 iJ 1 S 2 £1 

1: lia-'PL + ¥Pr) 1: ^(a-^PL + b^PR) 

2: i{b'*PL + a'*PR) 2: i{b'*PL + a'*PR) 

3: -ieqs{p-2kY 3: -ieqfY 



Figure 1: One- loop SUSY diagrams contributing to radiative, i^'j, decays. In (a), 

the scalar is a charged slepton and the fermion / is a neutralino. In (b), the scalar S is 
a sneutrino and the fermion / is a positively [negatively] charged chargino (g/ = 1). 

6a^ < 3 X 10^^. There is no experimental evidence of an nonzero EDM for the electron 
(de). The most stringent upper bound, obtained in ref. [44], is de < 1.6 x 10"^'' e cm at 
90% CL. Likewise, there is no experimental evidence for radiative flavor-changing charged 
lepton decays. The 90% CL upper limits to the branching ratios for the muon and tau- 
lepton radiative decays are given by: BR(yU 67) < 1.2 x 10^^^, BR(r — > 67) < 1.1 x 10^^ 
and BR(r ^ ^j) < 6.8 x 10"^ [11]. 



4.1 Supersymmetric corrections to the lepton-photon vertex 

The amplitudes for the processes of interest are obtained by evaluating triangle diagrams 
that contribute to the one-loop correction to the lepton-photon i'^i^'j vertex. Supersym- 
metric corrections to this vertex arise from the two topologies of diagrams depicted in fig. [H 
The corresponding Feynman rules required for the vertices are given in eqs. ( 1C.3I) and ( ]C.4I) 



of Appendix C The anomalous magnetic moment and electric dipole moment (EDM) of 
the leptons and the lepton flavor violating decays — >^ £^7 are derived from the following 
terms of an effective Hamiltonian: 

= e {Ci^ F a^""" Pd-^ + Cli' F a^"^ PrF) F^^ , (4.1) 

which can be extracted from the computation of the effective one-loop i^i'^j vertex. 

The computation of the Wilson coefficients Cl, Cr is straightforward. After calculating 
the contributions of diagrams (a) and (b) of fig. [T] and expanding in momenta of external 
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particles, we find for their total Wilson coefficients 

Cl' = C\'j( + m,.Cf/ + m,.Cf/ , (4.2) 

where the index i labels the contribution of diagrams i = a,b and the m^/ (/ = 1, 2, 3) are 
the lepton masses. For diagram (a) we obtain, 

Ql' = ^^,qsa'*a'C2s{ms, rrif) , = -l-^^q^hH' C,,{ms, mj) , (4.3) 

and for the diagram (b), 

C\l = j^<lfb'''a'mfCii{mf,ms) , Cl',^ = j^qfa'*b'mfCu{mf,ms) , 

C'd' = ^^g/a^*a^C23(m;, m^) , Cl^' = ^^^qjbH'C^simf, nis) , (4.4) 

where mj and nis are the masses of the fermion / and scalar S, respectively, and all other 
parameters are defined in fig. [H The loop integrals appearing in eqs. (14. 3 p and (14. 4p are: 

x'^ + y^ 2xV , y 

^^^^^'^^ = " 2(x^-y^)^ -(^^^^"gx' 

= - I2{x- -y- f + (^^^ X • ^'-'^ 

The full Wilson coefficients and Cr are obtained by summing over all relevant triangle 
diagrams in the model. In our case just two of them contribute: diagram (a) with charged 
slepton and neutralino exchange and diagram (b) with sneutrino and chargino exchange. 



4.2 {g - 2)^ and the electron EDM 

The formalism described above leads easily to expressions for the EDM of the electron and 
for the muon magnetic moment anomaly [g^ — 2)/2. For both processes I = J, so that the 
fiavor- diagonal piece of the effective Hamiltonian is given by 



^Ff"" , (4.6) 
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where we used the relation C/^ = C(l*. By matching to the standard form [45,46]o 

^ = --^ajFa,j'F^^ + '^Fa,,^,i'F^^ , (4.7) 
Arriij 2 

where aj = {gj — 2)/2 is the magnetic moment anomaly and d^j is the EDM of the lepton, 
one can extract the expressions for the electron EDM, df.-, and for — 2, 

4 = -2e\mC\l, (4.8) 
a, = -Am, [ReCll + m,{Cll + Cll)\ . (4.9) 

In principle, both quantities can be used to set bounds on parameters such as M, m^, 
m\ and Xj, that govern the heavy sneutrino sector. However, the one-loop contribution 
to the CW from fig. [TJ^b), which is sensitive to the sneutrino sector, is real if the chargino 
parameters /i and M2 are real. Hence, the electron EDM measurement does not yield 
any constraints on sneutrino parameters at one loop. However, there can be sensitivity 
due to potentially large two- loop corrections; for further details see Ref. [29]. Similarly, 
the neutrino magnetic and/or electric dipole moments^ are also insensitive to the heavy 
sneutrino sector at one-loop, since there is no possibility of attaching the photon to a 



one- loop graph that involves the sneutrino- neutrino-neutralino vertex (see Appendix C). 

The amplitudes displayed in fig. [T] can give sizable contributions to the anomalous 
magnetic moment of the muon. These contributions are flavor diagonal and are sensitive 
mostly to the overall mass scale of the sleptons, gauginos and light sneutrinos — i.e. to the 
diagonal entries of corresponding mass matrices. Thus, the measurement of a, can be used 
to set lower bound on these SUSY masses. Assuming that the discrepancy between the 
experimentally observed muon anomalous magnetic moment and the theoretical prediction 
of the Standard Model, 5a, ^ 3 x 10~^, is due to new physics effects arising from the 
diagrams of fig. [H one can deduce lower bounds on the magnitude of slepton squared-mass 
parameter as a function of M2 and tan j3. Examples of such bounds are listed in Table [21 

Note that potential contributions to Mf^^j [cf. f l3.39l) ] from the terms containing the Dirac 
mass niD are suppressed by a quantity of (9(m^M^^). As we will show in Section [^^31 this 
ratio can be at most of the order of 10~^, otherwise the Dirac mass term would generate 
unacceptably large contributions to rare i'j decays. Thus, the muon anomalous 

magnetic moment can be effectively used to set a lower bound on the diagonal 22 element 



""^^In eq. (14. 7|) . the unit of electric charge e is taken positive, so that the electron charge is — e (which also 
coincides with the convention adopted by refs. [45] and [46]). Eq. (|4.7p is consistent with the corresponding 
effective Lagrangian of ref. [45] , by noting that Commins et al. define the anomalous magnetic moment of 
the electron to be k = — Oe (J.D. Jackson, private communication). 

^^Note that for Majorana particles only transition dipole moments can be nonzero. 
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Ms = 100 M2 = 200 M2 = 300 



tan p 
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\min 
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r 
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iiU 




10 


300 


270 


210 


15 


420 


420 


370 


20 


530 


570 


530 


25 


650 


740 


700 



Table 2: Lower bounds on the square root of {m\)22 from the measurement of a^. All 
masses are in GeV. 

of the soft slepton squared-mass matrix m\ and on the gaugino mass parameter M2, as 
specified in Table |2J The dependence on and /i is significantly weaker. 

4.3 Radiative charged lepton decay: i"^ — ^ ^^7 

The i'-' i^'j decay width is given by 

r{i'^i'i) = "^{\cir + n'n . (4.10) 

The corresponding branching ratio is obtained by dividing the result of eq. (14.101) by the tree 
level decay width, r{i'^ C.^v^v^) = m^jG'jp /192tt^ (where we ignore VT-propagator effects 
and a very small correction due to the nonzero mass of the light final state charged lepton) . 
In particular, the branching ratios for the experimentally interesting decays /i — 67 and 
r ^ yU7 are given by: 

BR(/x ^ 67) 

and 

BR(r fi-f) 

At leading one-loop order, fig. [Tt^a) yields an amplitude that is proportional to the 
off-diagonal terms of the slepton soft mass matrix m|^, and thus not relevant for setting 
bounds on heavy sneutrino parameterJ^. The amplitude corresponding to fig. [11(b) depends 

^^Of course this diagram is relevant when yi,-dependent corrections to nij^ entries are generated by the 
renormalization group evolution of parameters. This effect has been studied extensively in the literature 
(see e.g., ref. [47]), and we will not repeat this discussion here. 



'iOii e /|^i2|2 I |^12|2\ 



(4.11) 



487rV 



(4.12) 
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directly on the lepton flavor conserving part of the light sneutrino mass matrix, Mf^^j. This 



can be verified by using the Feynman rules collected in the [Appendix C| and employing the 
mass insertion approximation (MIA) expansion; for more details see e.g. ref. [48]. Assume 
(at least formally) that sneutrinos are closely degenerate in mass, 



(4.13) 



and then expand the functions Cjf or Cj^ [denoted generically in eq. f l4.14p by /], which 
depend on the squared-massed m| , up to the first order. This results in 



df 



dm% 



m, 



df 



df 



dm% 



(4.14) 



where there is an implicit sum over k. The advantage of this procedure is that it allows one 
to perform the sum over the sneutrino flavor index k in evaluating eqs. (14.111) and (14.121) . 
For example, the neutrino squared-masses always appear multiplied by a pair of sneutrino 



mixing matrices (due to the form of the sneutrino couplings given in Appendix C). Using 



the inverse of eq. (13.751) . one obtains Zl^Zi'^m'^^ = {M.y^j ■ 

It is possible to relax the assumption of approximately degenerate sneutrino masses. 
In particular, it can be shown diagrammatically that it is better to use appropriate ratios 
in place of the derivatives of eq. (I4.14p in the MIA expansion. Thus, for J > I (corre- 
sponding to the decay of a heavier lepton i"^ into a lighter lepton i^) and neglecting terms 
proportional to the lighter lepton mass, one arrives at the simple result: 



C 



ij 







c 



IJ 



bIJ 



R — ^IR + ^l-' ^AL 
2 



(4vr)2 2.2^ 



[Mlc] 



72 



\ Am^ J .J J cos (3 Mw 



(4.15) 



where the Z± are the chargino mixing matrices defined in ref. [7] 

( Cjj (m^+ , m^i ) - Cjj (m^+ , j ) 



/ AC,, 



2 2 

m-j — m-j 



klj 



for / ^ J, 
for / = J . 



(4.16) 



and are the three "CP-averaged" sneutrino masses, given by the positive square roots 
of the eigenvalues of M^^ [cf. eqs. (K^ and (1338)) ]. 
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tan (3 



10 



20 



Mlc 



... > 2702 <312 
> 2702 



> 5702 < §2 < 452 \ 
... > 5702 <1502 
> 5702 



Table 3: Bounds on the structure of the matrix elements of for M2 = ^ = 200 GeV. 
All masses in the Table are given in GeV. 

Clearly, our approximate expression for C^^ given by eq. (14.151) . which enters the decay 
rates in eq. (14.101) . is proportional to the lepton number conserving squared- mass matrix, 
M|^, defined in eq. (13.391) . Even in the case where m\ is diagonal, contributions to radiative 
lepton decays arise from the off-diagonal elements of governed by the general form 
of the matrices mo and [cf. the third term in eq. (13.391) ]. Notice that the flavor 
dependence disappears completely in the limit of diagonal m\ and = in which case 
is diagonal. 

The effect of the seesaw contribution to the lepton number conserving part of the sneu- 
trino squared-mass matrix, M^^, has not been previously noticed in the literature. This 
yields an extra contribution to the decay branching ratios BR(£'^ — >■ £"^7). Consequently, 
for a fixed set of chargino sector parameters (/i, M2 and tan (3) and soft slepton squared- 
mass matrix {mj^), the experimental bounds on the radiative lepton branching ratios can 
be used [via eqs. (14.111) . (I4.12p and (I4.15p ] to determine upper limits on the off-diagonal 
matrix elements of M|^. Examples of such bounds for M2 = /i = 200 GeV and two sets of 
tan/5 and rn™*"^ (previously exhibited in Table [2]) are shown in Table [31 In obtaining these 
bounds, we assumed that m| is diagonal so that fig. [U^a) does not contribute to the decay 
amplitudeo We then varied the matrix elements of until the constraints from mea- 
surements were violated. Moreover, we incorporated the full numerical one loop calculation 
for — i> £^7, presented in Section Wl] rather than the approximate expressions given, e.g., 
in eq. (I4.15p . Notice that there exist lower bounds for the diagonal elements of M^^ from 
{g — 2)^, but upper bounds for the off-diagonal elements of from BR(£'^ — > + 7). 

The results of Table [3] illustrate that the bounds on the square roots of the off-diagonal 
elements of are at least 10 — 100 times smaller than the square roots of the diagonal 



-"^^Non- vanishing off-diagonal elements of mj^ should in most cases tighten the bounds on M^^,, barring 
accidental cancellations between the amplitudes obtained from fig. [ija) and (b) . 
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elements. It is convenient to rewrite eq. (13.391) in the following form: 

= ml + ^Ml COS2/3 + m}jM-^m%{t + M-'^m%)-^M-^ml + 0{v^M-'^) 

= ml + \Ml cos 2/? + m*j^M-^mljM-^m^j^ + Oiv^M-"^) + 0{v'^m%M-*) , (4.17) 

where we have expanded out the quantity (1 + M~'^m\i)~^ under the assumption that 
II M~^m^ II < 1 (to be justified shortly). Eq. (I4.17P implies that the off-diagonal elements 
of are roughly of order m^^m^/M^ (barring any accidental cancellations). If we 

assume that m/j is of order the electroweak scale, then the bounds on the off-diagonal 
elements given in Table [3] imply that 

-g^OdO-). (4.18) 

with the strongest bound given by /i ^ 67 decay. This result suggests that ||m^||^/^ cannot 
be larger than about 10% of the Majorana mass scale M. Hence, + m\ ~ and for 
the estimates of the magnitude of the entries of the lepton number violating mass matrix 
M|y in the next section we henceforth set m\ = 0. 



5 Neutrino masses and the lepton number violating 
parameters 

In this section we examine the constraints on the lepton number violating sneutrino squared- 
mass matrix M|y from our knowledge of the physical (light) neutrino masses and mixing 
angles. 



5.1 One- loop contributions to neutrino masses 

The effective operator that describes the light neutrino mass matrix is given by: 

- ^^^^ = iM,>,V + H.c. (5.1) 

Note that t'/t'/ is a AL = 2 operator, since it changes lepton number by two units. In 
Section \3A] we evaluated the tree-level contribution to M^^ [cf. eq. (13. 8p ]. However, one- 
loop contributions to the light neutrino mass matrix can be significant, and in some cases 
these can be as or more important than the tree-level contribution [18,28]. The dominant 
one-loop graph involves a loop containing neutralinos and light sneutrinos, as shown in 
fig. [2]^a). Due to the presence of the lepton number- violating sneutrino squared- mass matrix 
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(a) 




vv 



if °, AO, G 




Figure 2: One-loop corrections to light neutrino masses, (a) The loop consisting of light 



sneutrinos {Sk, k = 1 . . .6) and neutralinos {xl 



1 ... 4) is the dominant contribution. 



(b) The loop consisting of a neutral Higgs (or Goldstone) boson and a heavy neutrino 
contributes a relative correction to the light neutrino mass of at most a few percent. The 
contributions of the corresponding graphs (not shown) in which the light sneutrinos in (a) 
are replaced by heavy sneutrinos and the heavy neutrinos in (b) are replaced by light 
neutrinos are suppressed by an additional powers of 0{vM~^) as explained in Appendix D 



M|y, which violates lepton number by two units, fig. [2|^a) can contribute significantly to 
the light neutrino mass matrix. Other one-loop contributions shown in fig. [2](b), yield 
corrections to the light neutrino mass matrix of at most a few percent, and thus can be 
neglected. 

In order to establish the results just quoted, we begin by reviewing the relevant in- 
teractions that govern the one-loop contributions to the light neutrino masses. The light 
neutrino couplings arise from eq. (12.81) and the supersymmetric sneutrino-neutrino-neutral 
gaugino interactions. After isolating the interaction terms containing one neutrino field, 
one arrives at 



r2,,l r,J* 



-j={g2W^ - grB)vivi* + Y{.c. , (5.2) 



where and B are the SU(2) and U(l) neutral (two-component) gaugino fields, and 
g2 and gi are the corresponding gauge couplings. Using eqs. (13.41) and (13. 6p . it follows that 
z/^ ~ z/£ + m^M^^z/^ and z/£ ~ z/^ — M'^m^vi. Likewise, it follows from eqs. (13.331) and 
(13:371) that 



^1 ~* 



h 5 



vl - (M^ + mliY^Mm^^vi 



(5.3) 
(5.4) 
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Thus, the effective interaction involving (at least) one light neutrino field is given by: 
X, ^ -Yj'l^mi^y,* + vl^'^^Hl - imnM-Y' [m^f^l* + ^l^f Hi) 

+^{g2W^ - giB) [uy,* + nioMiM^ + m^*)-V/i>g + H.c. (5.5) 
v2 

In order to perform the explicit loop computations, it is convenient to rewrite eq. (15. 5p 
in terms of mass eigenstate fields. The Higgs field if| is expressed as [49] : 

Hi = ^[v2 + h° cos a + H^sma + i(cos pA^ + sin pG'^)] , (5.6) 

in terms of the CP-even Higgs fields hf^ and H'^ (where m/jo < tuho), the CP-odd Higgs field 
and the Goldstone field G", where tan/3 = V2/V1 and a is the CP-even Higgs mixing 
angle. We also define two-component mass-eigenstate neutralino fields (j = 1,...,4) 
following ref. [7] by 

tjji = Z%k] , where V'i = (-^^ , -iW'' , R\ , Rl) , (5.7) 

and Zjv is a unitary matrix that governs the Takagi-diagonalization of the complex sym- 
metric 4x4 neutralino mass matrix, M^o via Zj^M^oZj^j = diag(M^o , . . . , M^o). 

Before presenting the explicit computations, let us first estimate the order of magnitude 
of the loop-contributions to the neutrino mass due to the loop graphs of fig. [2](a) and (b), 
and the corresponding graphs (not shown) in which the light sneutrinos [heavy neutrinos] in 
graph (a) [(b)] are replaced by heavy sneutrinos [light neutrinos]. This analysis is presented 
Appendix D| — the results obtained there imply that the graphs of fig. [21(a) and (b) both 



m 



yield contributions to the one-loop light neutrino mass matrix of order the tree-level light 
neutrino masses, multiplied by the appropriate vertex couplings and a typical loop factor. 
Other one-loop contributions not shown in fig. [2] are suppressed by additional powers of 
(9(fM~^) and are utterly negligible. 

We begin with an examination of the loop amplitude of fig. MJo), which is governed by 
the light neutrino-heavy neutrino-Higgs interaction term of eq. (15.51) . The internal heavy 
neutrino line is marked with an x to indicate the lepton-number violating propagator 
proportional to its (diagonal) mass M6^^. Summing over all the internal neutral Higgs and 
Goldstone states, the leading 0{M) term vanishes, leaving a subleading term of 0{v'^M~^), 
which is the magnitude of the light neutrino mass. We find that fig. MJo) yields a leading 
contribution to the light neutrino mass that is proportional to the tree-level light neutrino 
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mass matrix [cf. eq. (13. 8p ]: 

where M and Wih denote average heavy neutrino and Higgs boson masses. This correction 
turns out to be of the order of at most few percent. Additional corrections can also arise 
that modify the flavor structure of M^,^, but these are not logarithmically enhanced and 
are thus even smaller. 

Hence, the possibility of a significant one-loop contribution to the light neutrino mass 
matrix can only arise from fig. [2](a), which is governed by the light sneutrino-neutrino- 
gaugino interaction term of eq. (15.51) . In the following, we examine the corresponding loop 
graph in which the external light neutrino fields are mass eig enstates (1//)?^^*^ [cf. eq. (ESI)]. 
Using four- component spinor methods, the amplitude for this graph (with incoming four- 
momentum p) will be denoted by 

- + 1:'s')Pl + m'v'* + *)P«] , (5.9) 

where the generic self energies T,yg{p^) of the Majorana neutrino must be symmetric in its 
indices /, J . To evaluate this graph, we express the neutrino-sneutrino-gaugino interaction 
Lagrangian in terms of the four-component self-conjugate Majorana neutrino fields and 



the Majorana neutralino fields [cf. Appendix A 
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= -1 {g^Z^^ - giZ]^){Zl' - iZ'i^'^'Mi^s ^P^mSu + H.c. , (5.10) 



where the neutralino mixing matrix Zn is defined in eq. (15.71) . The resulting DR- 
renormalized neutrino mass matrix at one-loop order is given by: 

(M(;-i°°P))" = + Re [E^/«) + lm.„S(/(m^J + lm.„S^^^«)] , (5.II) 

where the loop diagrams are regularized by dimensional reduction and the tree level di- 
agonal mass, m^^j, is defined at the renormalization scale hji. In addition, m^^, is some 
average neutrino mass scale, which to a very good approximation can be taken to be zero 
in the explicit loop calculations presented below. 

In order to determine the masses of the light neutrinos at one-loop accuracy, it is usually 
sufficient to calculate the diagonal matrix elements of the self energies (i.e., by setting I = J 
in eq. (15.111) ). assuming that the tree- level neutrino masses are non-degenerate. However, 



^^More explicitly, the non-zero components of Plv\[ are the two-component neutrino fields (t'|)P'^^'^, and 
the non-zero components of -PlX° are the two-component neutralino fields introduced in eq. (|5.7p . 
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in some cases S^'^/ can be numerically large for I ^ J . If the latter holds, then one 
must re-diagonalize the neutrino mass matrix, (M°^'^^~'°°p)'^"', in order to obtain the loop- 
corrected physical neutrino masses and corresponding mixing matrix [/mns (more details 
of a similar procedure in the context of i?-parity violating models can be found, e.g., in 
refs. [50] and [25]). 

An explicit calculation of the diagram shown in fig. [2]J^a), in the limit m^^ —>■ 0, yields 

V7J_~"M y2i ^ ^li\2f^Lk ■-;riL+3)k\ ( ^Mk ■ ^iM+3)k.-.j-LI jjMJ td(^ ^ \ 



4(47r) 



(5.12) 



V^/J_ ^ 1^ y2i ^^ry\i\2(rzLk ■ML+3)k. , ^Mk , ■'p-{M+3)k.-.^LI ttMJ*t^(^ ^ 



(5.13) 



with an implicit sum over repeated indices, where m^o and nis^ are the neutralino and 
sneutrino masses, respectively, and So, -Bi are the standard 2-point loop- integrals [51] 
evaluated at = 0, 

B,{x,y) = A-log^ + l-4^1og-, (5.14) 

\ 1 A 1 1 3 y"^ ( 1\, X , 

Bi{x,y) = --A + - log — - - - — ^ r- + r-^ - - log - , (5.15) 

2 2 yUjj 4 2[x^ — y^) \[x^ — y^y 2 J y 

with A = 2/ (4 — (i) — 7 -|- In 47r set to A = in the minimal subtraction renormalization 
scheme. Note that S5 is finite, i.e. in the sum over k the dependence on A and fiR cancels 
exactly due to the orthogonality of Z. Likewise, S(/ is finite for / 7^ J, which is easily 
verified after using the orthogonahty of Z and the unitarity of Umns- This is to be expected 
since in the mass basis there are (by definition) no tree-level off-diagonal neutrino mass 
matrix elements. In contrast, T^y^ is divergent, and after minimal subtraction it is here 
that the /ir dependence resides. 

We now examine the relative magnitudes of the various contributions in eq. (15. lip 
to the loop-corrected neutrino mass. First, we observe that Sy [given by eq. f l5.13p ] is 
dimensionless and has a magnitude of the order of a typical electroweak correction (this has 
been numerically confirmed). Thus, the one loop contribution of the terms proportional to 
the minimally subtracted Sy in eq. (15. lip is at most a few percent of the tree-level neutrino 
mass. Given the current experimental accuracy of neutrino data, this latter correction can 
be neglected, as it does not provide any constraints on sneutrino parameters. Thus, we 
focus on S5 [given by eq. (15.121) ]. which can be simplified by employing the MIA expansion 
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described in Section 14.31 The end result is: 



5M, 



ij 



1— loop\/J 
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(5.16) 



iKM 



where in analogy to f l4.16p we define 



V Am2 J 



m?7 



klj 



dBo{m rrij^i) 



for / ^ J, 



for / = J . 



(5.17) 



and the CP-averaged sneutrino masses, m^/, are defined below eq. (14.161) . As expected, 
this contribution is finite and is explicitly lepton number violating, as it is proportional to 
the matrix Mly. Eq. (l5:T6|l is a generalization of eq. (7) of ref. [18] to the 3-fiavor seesaw 
model, ^ 

The results given in Section 15.11 can be used to estimate the bounds on the heavy 
sneutrino soft parameters m^, m^, Xy imposed by the current experimental measurements 
of neutrino masses and mixing. These bounds allow for a significant one-loop correction to 
the light neutrino mass matrix, SM^^ , which could even compete with the corresponding 
tree-level masses. Further details will be given in Sections 15.31 and 15.41 



5.2 Radiative generation of neutrino masses and mixing 

It is very tempting to explain the characteristics of the neutrino mass spectrum as a conse- 
quence of radiative corrections. The most economical possibility is one in which the pattern 
of neutrino masses is entirely radiatively generated by the loop corrections. However, in 
the supersymmetric seesaw model this is not possible. If one sets mj,^^ = (for all I) in 
eq. (I5.16p . then m^) = (or equivalently, = 0), in which case only the light sneutrino- 
neutrino-gaugino interaction of eq. (15. 5p survives. However, this interaction generates a 
one-loop neutrino mass that is proportional to M|y [cf. eq. (15.161) ]. which vanishes in the 
limit of m/) = 0. 

Here, we shall be less ambitious and investigate whether the hierarchy and/or the fiavor 
mixing of neutrinos can be generated entirely by loop effects. As we shown below, such a 
scenario seems to be possible. However, in order to obtain the correct values of the light 
neutrino mixing matrix elements, a fine-tuning of sneutrino parameters may be required. 



We correct here a typographical in eq. (7) of ref. [18] where {g2Zf} — giZ]^Y is incorrectly written as 
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To be more specific, consider the following scenario. At tree level we assume the Yukawa 
coupling matrix Y^, to be real, non-negative and flavor diagonal, i.e. Yj'^ = YjS^'^ (with 

> 0). Consequently, the tree level neutrino mass matrix [eq. (13. 8p ] is also real, non- 
negative and diagonal so that f/j^NS ~ Then, the one-loop correction to the neutrino 
mass matrix [eq. (I5.16P ] is proportional to: 

«^^-^E-x?(^i^^'-^2^^T(^) ■ (5-18) 

i—l \ / ilj 

If one assumes that the flavor splitting of the light sneutrino masses is small, then the 
ratio {ABq/ Am'^)-jj is approximately constant with the respect to the indices J, J, so that 
tt/j a is roughly constant. Therefore, the one-loop corrected neutrino mass matrix 
[eq. (15. lip ] can be written as 

^a-ioop) ^ -mijM-^mD + Re {aMly) . (5.19) 

Since we have assumed above that Y^ is diagonal, it follows that m£, = V2Yy/\/2 is also 
diagonal, in which case there is no need to distinguish between m^) and its transpose. For 
simplicity, we shall further assume that <^ M^. Then, using eq. (I3.40p for M|y, in 
which only the leading 0{vM~^) terms are kept [under the assumption that m\ ~ 0{vM) 
as suggested by eq. (I2.14p ]. we may express eq. (I5.19P in the following form: 

^a-ioop) ^ _[i _ Re(aX^)] moM-^mD [1 - Re(aXj)] 

- 2mD^Re{aml)^mD + Re(aX^) mDM"^mD Re(aXj) . (5.20) 

To achieve the correct hierarchy of neutrino masses and mixings, one possible strategy 
is to demand that the sum of the last two terms on the right hand side of eq. (I5.20p is 
negligible, in which case the first term yields the correct physical neutrino masses and the 
mixing matrix. Then, using eq. (12. 5p . we perform a Takagi-diagonalization to identify the 
physical (loop-corrected) neutrino masses and mixing matrix elements: 

- [1 - Re(aX,)] mnM-'mn [1 - Re(aXj)] = (f/^^)* mj^^ ([7^^)+ , (5.21) 

where m^y'^ is the (non-negative) diagonal physical neutrino mass matrix. One can solve 
eq. fl5.2ip analytically for Re(aX,^), which yields: 

Re(aX,) = 1 - i{uZin<!:n"^RM"^m-^' , (5.22) 

where i? is a complex orthogonal matrix, subject to the restriction that the right hand 
side of eq. fl5.22p is real. Thus, starting from any hierarchy of the tree-level diagonal, non- 
vanishing Yukawa couplings yj, the special choice of X^, given in eq. fl5.22p allows us to 
reproduce the correct neutrino mass hierarchy and the mixing matrix. 
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Clearly, the scenario just presented is not very realistic from the phenomenological point 
of view. To achieve the desired result, a specific form of the parameter, very close to 
perturbativity limit of Y,y and the charged slepton masses is required, as well as a rather 
precise cancellation between the last two terms of eq. (15.201) . Nevertheless, our example 
above provides an analytical existence proof for a radiative mixing scenario. In general, 
for given Y^, and M, many choices of sneutrino parameters leading to the correct pattern 
of neutrino masses and mixing at the one- loop level exist, but they need to be determined 
numerically. Presumably, all successful scenarios require a certain degree of fine-tuning, 
but perhaps some solutions would be deemed acceptable. 

5.3 Universal parameters at the scale M 

The magnitudes of the parameters A^, m\ and that govern the behavior of the 
heavy sneutrino sector are connected with the mechanism of supersymmetry breaking [cf. 
eq. ( 12. 6p ]. These parameters decouple at the scale M ^ Mz where the sneutrino superfield 
N decouples. If the scale M is close to the GUT scale then soft SUSY breaking parameters 
are restricted by GUT symmetry considerations. Further assumptions on the minimality 
of the Kahler potential in supergravity simplify our input parameters considerably, at the 
scale M ~ Mqut, 

Ay = AqY^ , ml = moM , = xM^ , (5.23) 

where Aq is a complex number, mo and x are real numbers, M is a diagonal 3x3 Majorana 
neutrino matrix [cf. eq. (12.51) ] and Y^, is the neutrino Yukawa coupling [cf. eq. (12. ip . 

Under the universality assumptions of eq. (15.231) . the matrices M^^ and M|y assume 
the following simple forms at the GUT scale: 

Mlc = < + ^Mlcos2/3 + -^m},ml, (5.24) 

Mly = ^(^o + /.*cot/5-^) , (5.25) 

where the light tree- level neutrino mass matrix Mj,^ is given in eq. (13.81) . As parameters 
"run" from the GUT scale to low energies, m| receives renormalization from other Yukawa 
and gauge interactions. In contrast, all the parameters associated with the superfield N 
are hardly affected since M ~ Mqut- Moreover, the neutrino mass matrix M^,^ and the 
superpotential parameter /i are both multiplicatively renormalized. Hence, just above the 
scale of low-energy supersymmetry breaking, the low-energy value of Mly is still given by 
eq. (I5.25p . with the parameters on the right-hand side defined at the low scale. At the 
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low-energy supersymmetry-breaking scale the DR running neutrino mass matrix M^^^^r) 
[or its diagonal form m^7(/iij)] receives finite threshold corrections from the neutralino- 
sneutrino loop in fig. [2](a). The one-loop correction to the neutrino mass matrix given in 
eq. fl5.16p is proportional to the diagonal tree-level neutrino mass matrixo Hence, the 
one-loop corrected neutrino masses assume the very simple and suggestive form 



(l-loop) 
^1,1 = ^u,l 



l + 2Re . " . lAo + fi*cotp 



1 + x) V 1 + X 



(5.26) 



where a is defined in eq. (15.181) and all parameters are now defined at the scale fiR = Mz- 
We next examine the light sneutrino mass difference. Since the results of Table [3] imply 
that is very close to diagonal form, it follows that — 1 (cf- discussion above 

eq. (13.481) ]. Combining the results of eqs. (I3.49p . (I3.52p and (I3.66p . we derive 



' Ao + /i cot fJ - 



1 + X V 1 + a; 



(5.27) 



which is identical to the one flavor case found in eq. (13.681) and in Ref. [18] if the neutrino 
mass matrix M,y^ is diagonal. In the more general case of non-diagonal M^^, the diagonal 
elements of the neutrino mass matrix do not coincide with the neutrino masses m,. , . Con- 
sequently, the quantity {^rriyjmu^) ^ exhibits non-trivial dependence on the flavor index /. 

To produce quantitative results, we need to initialize the neutrino Yukawa couplings in 
such a way that we always reproduce the "observed" MNS mixing matrix. Using eqs. (13. 8p 
and (I3.10p . it follows that 

mo = ^t/^NS «^^)'/' M'/' , (5.28) 

where R is an arbitrary complex orthogonal matrix [47], with three (complex) angles, 6'i 2,3. 
(As the sign of R is undetermined, one may choose det R = 1 without loss of generality.) In 
the plots that follow, we assume a hierarchical spectrum for the neutrinos, and all relevant 
input parameters are displayed in TableHl The value for m^, adopted in TablelHis consistent 
with a supersymmetric interpretation of the observed experimental excess for Sa^. 



^"'^ Indeed, assuming universal parameters at the GUT scale, and noting that x < 0(10^^) [cf. eq. (|4.18[) ]. 
it follows that M^Q ~ at the GUT scale, where mj^Q is one of the approximately degenerate 

eigenvalues of M^^. The positive square roots of the eigenvalues of Ml^, evaluated at the low-energy scale, 
are identified as the three CP-averaged light sneutrino masses. Although nij^ is no longer proportional to 
the identity matrix at low-energies, this latter effect is formally of higher order in the loop expansion of 
SM^^ [cf. eq. (|5.16p ]. Consequently, we can neglect the flavor splitting of the CP-averaged light sneutrino 
masses in the evaluation of the ratio (ABo/Am^)^;^ in which case this ratio is roughly constant with 
respect to the indices K and M as discussed below eq. ()5.18p . 
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Input Parameters 



Neutrino Sector SUSY Sector 



^phys 


10-14 


Ao 





^phys 









^phys 


/ 




350 


ei 


0.2+0.1 i 


tan P 


10 




0.3 


Mb 


95 


O3 


0.1 + 0.5i 




189 


Ml 


10^4 


X 


0.0 


M2 


2 X 10^^ 


rriL 


197 


Ms 


5 X 10^^ 


rriR 


135 



Table 4: If not otherwise indicated, the input parameters that govern the neutrino and 
SUSY sectors listed above have been employed in our numerical analysis. We take Amg^j = 
(8.0i[!:^) X 10-5 eV^ and Am^^^ = (2.45 ± 0.55) x lO'^ eV^ from Ref. [11]. The values for 
6'i,2,3 above are representative choices (as these angles are not fixed by the light neutrino 
data). All mass parameters in the above table are in GeV units. 

In fig. [3] we plot the ratios {Anij^Jm^Jj [upper panels] and (mi], / fn,y^) j [lower 
panels] as functions of the SUSY-breaking parameters mo [left panels] and Aq [right panels] . 
When varying niQ we set Aq = and when varying Aq we set mo = 0. Otherwise, our input 
parameters are as specified in Table |H In obtaining these results, we have incorporated the 
full one-loop contribution to the neutrino masses. In the two lower panel plots, the ratios 



m 



(l-loop) 



are nearly independent of the flavor /, and thus only one curve is shown. 
Our numerical results confirm our analytical approximate formulae of eqs. fl5.26p and fl5.27p 
and demonstrate that one must have mo < 10^ GeV {\Aq\ < 10^ GeV) to guarantee that 
the radiative corrections to neutrino masses are less than 80% of the tree level neutrino 
mass. In this case, the sneutrino mass difference is at most Amc^ < 300 Am atm — 15 eV. 

For completeness, we plot in fig. H] the results for — 2 anomaly and the branching 
ratios for the decays f'^ i^j in the case of universal parameters at the SUGRA scale. 
The results shown in fig. H] confirm our choices of a lower bound for mi [cf. Table |2] obtained 
in Section Wl] and an upper bound for x [cf. eq. (14.181) ] obtained in Section l473l 
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Figure 3: Predictions for the ratios {Arriyjmy^) j and {ttl'u^^'^^^'' / rrijj^) ^ for the three neutrino 
states (/ = 1, 2, 3) as functions of the soft SUSY-breaking parameters mo and Aq. When 
varying itiq [left panels] we set = and when varying Aq [right panels] we set mo = 0. 
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Figure 4: (a) In the left panel, the contribution to the muon anomalous magnetic moment 
from the diagrams in fig. [1] as a function of mi = mji is exhibited, (b) In the right panel, 
the prediction for BR.{i-' — * i^'j) is shown as a function of the parameter x = mj^/M'^. 
The upper [lower] curves correspond to r — > /i7 [r — 67], and the middle curve to /x — 67. 
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5.4 General case 



So far we have dealt with universal boundary conditions for the supersymmetric parameters. 
One can set general bounds for the lepton number violating matrix elements of M|y from 
eq. (15.161) and the "naturalness" assumption of drriy^ < rriy^. In the general case, appropriate 
bounds can be derived only numerically and depend on the particular form of the MNS 
matrix. Analytical estimates can be obtained using the following approach. Let us require 
that the one-loop corrections to the neutrino mass matrix do not significantly affect the 
physical neutrino masses and their mixing. Combining Eqs. fl3.10p and fl5.16p . one gets for 
any /, J: 



"^^i -D^\!., ^ nrU\2TTMI ttNJ //i^2 \ 1 (^Bq 



,2 



iMN 

(5.29) 



The structure of the ?7mns factors on both sides of eq. fl5.29p is identical, so roughly [barring 
possible cancellations between terms and the effects of truncating a potential imaginary 
parto of U^^^ {.^i^i)mn ^mnsIj condition above can be rewritten as: 



^ \aMN{Mly)^^^^\, (5.30) 

with aMN defined in eq. (15. 180 . 

Further estimates depend on the particular choice of the m^j (or Fj,) and M and on the 
neutralino sector parameters. For example, using the parameters specified in Table HI one 
has aMN ~ a ~ 4 X 10"^ GeV~\ so that 

I [Mly)^^ I < 2.5 X 10^ GeV | (M,J^,^ | . (5.31) 

Eq. (15.3 ip implies that in the general case one should expect the entries of the matrix 
M|y to be no more than 5 or 6 orders of magnitude larger then the typical scales in the 
effective neutrino mass matrix; i.e. of the order of a few MeV^. Bounds on M|y can be 
also translated into bounds on Xy and m\. From eq. (l3.4Up one can see that, barring fine 
tuning, we have approximate relations ~ M^^Xu or M|y ~ My^m\/M. Thus the 

rough estimates we made above suggest that both X^ and m\/M should be smaller than 
approximately 100 TeV. 



iMN 



^^If the Higgsino mixing parameter /i and the lepton trihnear couphng Ai are real (the case of complex 
(1 and Ai has been extensively discussed in the literature, see e.g. [52]) then there is no bound on the 
imaginary parts of the matrices and 
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Stronger bounds on the matrix elements of M|y can be obtained numerically after 
assuming some particular form of the MNS matrix. As an example, under the assumption 
of tri-bimaximal mixing of ref. [53] and the parameters given in Table HJ 



< 



( 



\ 



2 X 10-9 



2 X 10- 



10" 



GeV^ 



(5.32) 



where the dots indicate elements with similar bounds as the diagonal ones. The significant 
suppression of the lepton number violating matrix elements of M\y relative to the lepton 
number conserving matrix elements M\q ~ Oiv^^ is particularly noteworthy. 



6 Sneutrino Oscillations 

The theory behind sneutrino oscillations follows closely the very well known theory of 
oscillations in the neutral Kaon- meson system. The light sneutrino state [cf. eq. (15.31) ]. 
i/£ ~ ~ m^M(M^ + rnj^)~^u^ is to leading order in vM~^ the supersymmetric partner 
of left-handed neutrino ul, and therefore couples to the and Z gauge bosons. For the 
present discussion, it suffices to approximate: z7/ — ^l-, which we shall denote simply by vj 
in this Section. The vi can be produced, for example, in e^e~ annihilation via s-channel 
Z exchange: 

e+ + e-^Vi + V}. (6.1) 

When lepton number is conserved, the z// (z/|) possess a definite lepton number equal to 
— 1 (+1) and they are produced in definite flavor eigenstates / = 1, 2, 3. 

It is convenient to introduce a two-dimensional complex vector space spanned by a 
basis of vectors consisting of the sneutrinos states of a given flavor /, |z//) and Two 
important operators that act on this state are: 




where L is the lepton number operator and CP is the CP-operator in the \ basis. 

That is, and |z/|) are eigenstates of L: 

L\vi) = -\vi) , L\v*)=+\v}), (6.3) 
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and the charge-conjugate parity operator CP transforms particle states into antiparticle 
states: 



CP\uj) = \i)*) , CP\i)*) = \i)i). (6.4) 
The eigenstates of CP are given by 

= + , 1^}"^) = T^d^/) - i^D) , (6.5) 

with definite eigenvalues 

CP|I7f+)) = , CPlu'j-^) = -lu'j'^) . (6.6) 

The CP-even sneutrino state of flavor /, l^'j^''), and the CP-odd sneutrino state of flavor 
/, \uj ■*), are states of indefinite lepton number. Of course, these states are the real and 
imaginary parts of the sneutrino field of definite lepton number, 

uj = ^{^^^ + . (6.7) 

Inevitably, in a supersymmetric model with a mechanism that yields neutrino flavor 
oscillations, the sneutrino flavor states should oscillate as well. The sneutrino mass eigen- 
states, Sk, {k = 1, 2. ..6) are linear combinations of the CP eigenstates and for a three 
flavor system (J = 1,2,3) they are related by: 

|z7W) = Zf|S.) , \^-^) = zi'^''>'\S,) , (6.8) 

where the real orthogonal 6x6 matrix with Z^^ has been introduced in eq. fl3.75p . The 
l^fc) are states of definite CP unless the following CP- violating conditions hold: 

Zi'-^^^^ , Z^^+3)J _^ Q ^ /, J = 1, 2, 3 . (6.9) 

In the presence of complex parameters in the Lagrangian (whose phases cannot be absorbed 
by field redefinition), one expects the conditions specified in eq. (16. 9p to be satisfied (even 
in the case of a one-generation model). 

Let us initially focus our analysis on the CP-conserving one-generation model. Consider 
the time evolution of the sneutrino states. The time dependence of a sneutrino in the state 
Iz/*-^) is governed by a definite frequency uj± = E±/h where E± = [jp'c^ + m^c^)-*^/^. where 
m+ and m_ are the masses of l^'*-^-') and l^'*-"-') respectively. If these masses are large 
compared to momentum p then the corresponding energies are E± ~ m±c^ (in which case, 
uj± ~ m± in units where h = c = 1). In addition to the time-dependent phase, we must 
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also account for the fact that the sneutrinos decay exponentially (e.g. into a chargino and 
a lepton) with a lifetime of t± (for respectively). We exhibit this time dependence 
explicitly by writing 

—iuj^t- 



27r|^{-)) 



(6.10) 



where the z/^^ are time-independent state vectors, That is, starting at t = 0, the proba- 
bility for finding particle in the sneutrino state u^^^ is given by |(z/^+)|\E'(t))p = 6"*/"^+, as 
expected. 

The well known striking effects of the K-system (e.g., K-K mixing and regeneration) 
can also occur in the sneutrino system. For example, we demonstrate how sneutrinos 
states can turn to states If we start off with a sneutrino state that is ^^(0) = |z/) = 
-^{\^^^^) + ■*)) at t = 0, then it follows that at time t. 



Then, the probability amplitude that the sneutrino is in state |z/*) is 
P,^,,{t) = \{v*\^>{t))\' = 



2e 



cos[(co'+ — ct;_)t] 



(6.11) 



. (6.12) 



The quantum interference effects can only be seen if t 
{Am)t = 0(1). That is. 

Am 



r_ and (m+ — m_) t 



r,- 



0{1) 



(6.13) 



where F/j is an average decay rate for the sneutrino, and Am is the mass difference of the 
CP-even and CP-odd sneutrino states. Eq. f l6.12p describes the oscillations of sneutrinos 
into antisneutrinos, or equivalently the oscillation between states of definite CP quantum 
number. We shall call this phenomena CP-driven oscillations. 

Similarly, one may compute the probability that the initial state \i>) is in the state \i>) 
at time t. We find 

1 



P,^,(t) = |(z>|vl/(t))|^ 



-t/T+ 



+ e-*/^- + 2e cos[C 



Co" I 



(6.14) 



One can also easily verity that Pj:,*^^* = Py^i, and Py*^y = Pj;,^^*. However, the proba- 
bility Pu^u is proportional to the number of negatively charged leptons (A*";-) due to the 
decay z/ — ^ + while Pa^c* is proportional to the number of positively charged leptons 
(A*';+) due to the decay z/* — > /+ + Then the asymmetry. 



A, 
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(6.15) 



is proportional to the quantum interference term cos(Amt) in eqs. (16.121) and (16.141) . That 
is, the lepton charge asymmetry Ai oscillates in time and provides a possible method for 
experimentally determining the value of Am. 

The signal for sneutrino-antisneutrino oscillations can be interpreted as the observation 
of a sneutrino that decays into a final state with a "wrong-sign" charged lepton. The 
phenomenological implications of such wrong-sign charged lepton final states at future 
colliders have been explored recently in Ref. [54]. 

We now turn to the three-generation model (allowing for the possibility of CP-violation) 
and consider the additional possibility of flavor metamorphosis. We pose the following 
question: Given the state {uj) at time t = 0, what is the probability that the sneutrino at 
time t is in the state |z/}) or |z/j)? Following the arguments given above eq. (16. lip , we find 
that a sneutrino wave function involves with time according to 

|v[/,(t)) = -L (Zf + iZi'"-'^') e-'^''-^ \Sk) . (6.16) 

Hence, the probabilities to be in the state or |z/j) at time t are given by: 

1 ^ [ 1 1 1 
Puj^u^it) = P^;^^,(t) = - e"H^+^J COS [(cufc - cu,)t] X 

k,s=l 

Z-^ Z-^ Z-'' Z-" + 2,^^+"^)^ 2,^I+'i)k 2(J+3)s ^(/+3)s _ 2,1^ ^(^+3)8 ^(7+3)5 

+ Z-^ 2U+3)fc 2 Js ^{/+3)s _^ 2{ J+3)A: ^/fc ^(^+3)5 ^Is _^ 2Z~^ ^(^+3)^ ^{7+3)8 ^Is\^ 

(6.17) 

1 ^ [ 1 1 1 

k,s=l 

( T-^fc T-^fc T-Js cr/s I ^^(^+3)^ -r{^+3)fc MJ+'i)s ^{/+3)s 07 -r(<^+3)5 o-(/+3)s 

+ Z^'' 2^^+^)'' ^{^+3)8 _j_ J+3)fc2/fc J+3)s _ 22i^ 2U+3)fc ^{^+3)5 ^/s j 

(6.18) 

Note that the probabilities in eqs. (I6.17P and (I6.18P are unchanged under the interchange 
of flavor indices I and J, respectively. The three-generation model possesses both flavor 
and CP-driven oscillations. 

In the supersymmetric seesaw model, neutrino mixing and masses are governed by a 
variety of parameters that contribute to the tree-level and one-loop neutrino mass matrix 
(cf. Section [5.2p . Some of these parameters also are relevant for determining the struc- 
ture of the real orthogonal sneutrino mixing matrix Z^^ ^ which controls the properties of 
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the sneutrino mixing as shown above. Consequently, the bounds on the model parame- 
ters discussed in Sections H] and [5] can be used to significantly constrain the general form 
of eqs. (EUD and (lOSD . 

The mass splittings among sneutrinos of different fiavors is typically much larger than 
the sneutrino-antisneutrino mass splitting between sneutrino states of a given flavor. In 
particular, due to the renormalization group evolution of parameters, Amfj is generally 
larger than few GeV^, even in the case of universality assumptions at the high scale, whereas 
sneutrino-antisneutrino mass splittings are typically of order the light neutrino masses. The 
observability of oscillations depends on the ratio Am/F [cf. eq. f l6.13p ]. Because the total 
decay width, F, is universal for a given sneutrino, whereas the scales of the corresponding 
mass splittings are so different, it follows that Am/F ~ 0{1) can be satisfled only for one of 
the two oscillation phenomena. That is, at most one oscillation phenomenon, either flavor 
oscillations or CP-driven oscillations, can be observed. 

Consider flrst the CP-driven oscillations. These oscillations can be observed if the 
lifetime of the sneutrinos is sufficiently long (the appropriate numerical requirements are 
given later in this section). In this case, fiavor-driven oscillations are much faster and 
have a very short "baseline" , so these oscillations are unobservable in collider experiments. 
Therefore, one can take a time average over fiavor-changing terms in the sums in Eqs. (I6.17P 
and (16.181) . setting them effectively to zero, and retain only those terms where the mass 
splitting is CP-driven and not flavor-driven (i.e. keep only those terms with s = k or 
s = k + 3). Now, the sum over s can be performed, and eqs. (I6.17P and (16.181) simplify to: 

3 

K=l 



K=l 

3 _ [ 1 1 ' 

- 2 ^ e'l^^^i cos[AA^t] Re{X^^X^^Y^^Y^^) , (6.19) 

K=l 
K=l 

+ 2 e'te"""^] cos[A^t] Re (X^^X^^y^^r-^^) , (6.20) 



K=l 



where A/< = luk — ^^k+s and we have used eq. (13.791) to express the 6x6 matrices in 
terms of the 3x3 matrices X and Y. 

Eqs. (16.191) and (16.201) are easily interpreted. For "long baseline" oscillations, one needs 
flrst to project flavor / onto some K (via the X^^,Y^^ factors), then the CP-driven 
oscillation takes place between the would-be sneutrino-antisneutrino states Sk and Sk+3, 
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and finally the result is projected back onto flavor J. 

Further simplification is possible if we exploit the bounds on the parameters due to the 
i-^ ^ decays obtained in Section H73l to conclude that the matrix Ml^ is very close to 
diagonal form. In this case, the matrix Qo that diagonalizes [cf. eq. (13.471) ] is close to 
the identity matrix. Moreover, the matrix elements of R [cf. eq. fl3.59p ] are suppressed by 
the ratio of Arrij^/mi^, and are therefore negligible. It then follows that X ^ F ^ T/v^, 
where T = diag(e~^'^i/2^ e-''^^^'^ , e'"'^-'/^) and (pj ~ axg{Mly)jj [cf. eq. fl330D ]. If we 
consider flavor conserving (i.e. I = J) sneutrino-antisneutrino oscillations, then there is 



one large contribution in eq. (16.191) in the sum over K for I = K, whereas the contributions 
oi I ^ K are strongly suppressed by the squares of mixing angles. Therefore, the dominant 
contribution to the probability for sneutrino-antisneutrino oscillations is given by: 



Pi 



1 

4 



+ e 



-t/TI_ 



2e 



i_ 



1_ 

-2ti_ 



cos(A/t) cos(207 



(6.21) 



which coincides exactly with the formula obtained previously for the one generation case 
[cf. eq. (16.121) ] in the CP-conserving limit (where Mf^y is a real matrix so that cos 20/ = 1). 
Similarly, for Pc^^c^, one reproduces eq. (I6.14p in the same limiting case. 

To complete the analysis of the sneutrino oscillation formulae, we must compute the to- 
tal sneutrino decay width, = T{Sk anything) = . Supposing that the neutralino 
is the lightest supersymmetric particle (LSP), the sneutrino decay width is the sum of the 
partial widths of the following two kinematically available decay chains 

3/2 



2 ^S, 

y2 



T{Sk 




111 2 



z 



Ik\2 



+ 12. 



a+3)fc|2 



(6.22) 



% 647r 



^2i„ 1 2 



E 

j=i 



- tZ. 



(J+3)fc 



JI 

MNS 



(6.23) 



In deriving the formulae above, we have used the Feynman Rules eqs. (IC.ip and flC.4p 
from Appendix C and have taken the lepton masses to zero. Eqs. (16.221) and (16.231) agree 
with Ref. [18] in the limit Umns = Zc = 1. Writing Z^ in terms of X and Y [cf. eq. (13.790 ], 
it easily follows that the decay rates of the sneutrinos Sk with = 1, 2, 3 [fc = 4, 5, 6] 



depend on X [Y] alone. Since X and Y differ only by the "small" R matrix [cf. eq. (I3.60p ]. 
it follows that r/_^ ~ r/_, which can be used to further simplify the expression given by 
eq. 



23r(S'fc ^ i"^^ + xf) indicates the sum of the sneutrino partial widths to the lepton-chargino and its 
charge-conjugated final states. 
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The total sneutrino decay width is given by: 

3 2 3 4 

7=1 i=l 1=1 i=l 

3/2 4 / ™2 \ 3/2 



2^ 

^2 327r 



(6.24) 



where the summation over the lepton indices can be performed in the hmit of vanishing 
lepton masses, with the use of the orthogonahty [unitarity] relations for the matrices 

[f^MNs]- 

How can one observe sneutrino CP-oscillations? Consider the following scenario: sup- 
pose that the LHC finds sneutrinos with masses that are accessible at a future Inter- 
national Linear Collider (ILC). Then, at the ILC, the sneutrinos are produced through 
the annihilation process of eq. (16.11) . and subsequently decay into [leptons + charginos] 
and [neutrinos + neutralinos] following the decay widths given by eqs. (I6.22p and (16.231) . 
respectively. Sneutrino CP-oscillations will then be observed only if the asymmetry Ai 
defined in eq. (16.151) . is appreciable, i.e., Ai ~ C'(l), which can be realized if both Am^ 
is small (providing a long enough oscillation base) and the sneutrino decay rate is suffi- 
ciently slow such that Arrik/Tk ~ C>{1). This scenario is impossible if the sneutrinos are 
sufficiently heavy compared to the neutralinos and/or charginos, in which case (neglecting 
the phase space suppression in eq. (I6.24p and performing the summation over the chargino 
and neutralino indices) the sneutrino decay rate is approximately given by: 
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(6.25) 



The expression above depends only on the sneutrino mass and cannot be suppressed by a 
particular choice of mixing angles of the Zj^, Z+ or Zn matrices. Thus, using the results of 
Section [5l one can check that the ratio Aruk/Tk is always much too small for the sneutrino 
oscillations to be observed. As an example, in the case of universal parameters discussed 
in Section [531 for the lightest sneutrino and mo, |^o| ^ 10^ GeV we obtain 

^ < 2.7 X 10-6 , (6.26) 
J- s 

which is very far from the value 0{1) required for the observability of sneutrino oscillations. 



48 



In the case of 2-body decays, the decay width Tk can be only suppressed by choosing 
an appropriate hierarchy of particle masses. Most of the decay channels in Eqs. f l6.22p and 
f l6.23p would have to be closed kinematically, with the open channels strongly suppressed 
either by the very small phase space factors (which requires rather unnatural degeneracy 
between sneutrino and neutralino or chargino masses), or by sufficiently small mixing angles 
for the relevant channel. An alternative possibility is one where the sneutrinos are lighter 
then all charginos and neutralinos, so that all 2-body decay channels are closed, but heavier 
than some charged slepton. In this case, z/ — >• i'^W'^, and assuming that the W is produced 
off-shell the end result is a 3-body decays that can produce an observable charged lepton. 
Three-body phase space significantly suppresses the sneutrino decay rate (relative to the 
two-body decay rates discussed above), and can yield observable sneutrino-antisneutrino 
oscillations, as shown in ref. [18]. However in such a scenario, either the charged slepton is 
the LSP, which is strongly disfavored by astrophysical data, or the charged slepton decays to 
some new lighter supersymmetric particle, which requires extending the model beyond the 
seesaw-extended MSSM considered in this paper [55]. As we have shown, the oscillations 
in the three-generation case does not differ much from the one-generation case, where the 
fiavor indices are summed over [cf. eqs. fl6.2ip and (16.240 ]. Thus, the results of ref. [18] 
can also be used without significant changes in the three-generation case discussed in this 
paper. 

Finally, we discuss the case of sneutrino flavor oscillations. These oscillations are de- 
scribed by eqs. f l6.17p and f l6.18p with indices I ^ J. For any choice of / 7^ J, both 
equations can be significantly simplified using the bounds on the structure of sneutrino 
mixing matrices derived in Sections H] and [51 These bounds imply that the off-diagonal el- 
ements of matrices Q and R [defined in eqs. (13.491) and (I3.56p ] are small, which then imply 
[via eqs. (13.601) and (13.791) ] that the off-diagonal elements of the matrices X, Y and are 
likewise small. Thus, to a good approximation one can keep in eqs. (I6.17P and (16.180 only 
terms at most quadratic in the non-diagonal elements of Zj^. For example, in the sum of 
the first term of the product of four Zj^^s in eq. (I6.17p . it is sufficient to keep only terms 
with s. A; = J, J -|- 3, J, J + 3. Assuming that the lifetimes of all eigenstates are very similar 
(i.e., r ~ r,fc), all the dominant terms can be summed to give a simple final expression valid 
for / ^ J , 

Pu,^uj ~ e-^ jig^g^-^l' + \Q''Q''*\'' + 2Re {Q''Q''*Q''*Q'') cos Am.jtj , (6.27) 

where Amu = f^ui — f^vj ■ 

The analogous expression for the sneutrino-antisneutrino oscillation probability Pi>j^-_»/>* 
is bilinear in the matrix elements of R [cf. eq. (13.560 ]. The latter are at most of 
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O{10^mi,/v) < 10 ^ and thus lead to completely negligible sneutrino-antisneutrino transi- 
tion ratesl^ 

The form of eq. f l6.27p is explicitly invariant with respect to rephasing, Q^"^ Q^^e^'f'-' . 
Thus, without loss of generality, we may replace Q by Qo [cf- eq. (13.491) ] in eq. (16.271) . 
where the off-diagonal matrix elements of the unitary matrix Qq are given approximately 
by eq. (13.481) and the diagonal elements of Qq are fixed by unitarity. As Qq is close to the 
identity matrix, the following approximations are valid: Qg"^ ^ 1 and Qq^* — —Ql/ for 
I ^ J . In this approximation, eq. (I6.27P simplifies for I ^ J to: 

P,,-,,, ^ 2e-*/" [\Ql''\^ -Re{Ql-^f cos Ami jt] • (6.28) 

If one uses the approximate expression given in eq. (13.481) . Ql/ ~ (M|^)^'^/(m?^ — m?^), 
then eq. (16.281) yields the oscillation probabilities directly in terms of the sneutrino squared- 
mass matrix elements. As expected, the sneutrino flavor-transition depends on the flavor- 
conserving matrix M|^. 

Defining the oscillation length hy L = ct we can write 

Am/jt = 5.06 X Am/j (GeV) L(fm) . (6.29) 

As in neutrino oscillations, it is useful to define Am/jL = 2ttL/Lq where Lq is the charac- 
teristic length of the oscillation : 

If the sneutrino mass difference is of (9(1 GeV), the characteristic oscillation length is of 
order 1 fm. Of course, the characteristic length of oscillation must be smaller than or at 
most comparable to the decay length of the particle for oscillations to be observable. In 
the case of the sneutrino, the decay length is [using eq. (I6.25p ]: 

28 (fm) 
mc, (GeV) 

Hence, the condition Ly > Lq requires that 

Amjj - 1 



> 

rui^ 25 



(6.32) 



^''An accurate estimate of Poj^u) should also take into account similarly small effects produced by the 
admixture of the heavy sneutrino states in the definition of the vi, which were neglected in derivation 
of eqs. (|6.17p and (|6.18p . However, given the extremely small transition probabilities, we do not present 
the full analysis here. 
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Such a mass splitting between the sneutrino states of different flavors is sensible. Thus, 
the likelihood of observing flavor sneutrino oscillations at colliders depends primarily on 
the degree of suppression caused by the mixing angles in the matrix Q. It is instructive to 
input some representative numbers in eq. (16.271) . Thus, for Ami2 = 10 GeV, mc, = 270 GeV, 
tan /3 = 10 and taking into account the bounds of Table [3l we obtain for z/g oscillations 
at time t = r = [cf. eq. ( KM ]: 

Pp^^^, ^ 1.25 X 10-^ [1 - cos(Ami2r)] , (6.33) 

Thus, as a consequence of the bounds from neutrino masses and radiative flavor changing 
decays obtained in Sections H] and [3, we conclude that in the see-saw extended MSSM, 
sneutrino flavor oscillations are difficult to observe at colliders. 

If the bounds of Sections H] and [5] could be avoided, say with some cancellation mech- 
anism (which in the absence of such a mechanism would appear unnatural), then it may 
be possible to find regions of the supersymmetric parameter space where flavor oscillations 
are observable. Then, at the ILC, one can define a flavor asymmetry for the number of 
muons vs. electrons in the final state, analogous to eq. (I6.15p . A time- variation of this 
flavor asymmetry would indicate the presence of flavor oscillations. 

7 Conclusions 

In this paper, we have studied sneutrino mixing phenomena in the seesaw-extended MSSM, 
allowing for the full complexity of the three-generation model (which includes both flavor- 
changing and CP- violating effects). We have focused primarily on the soft-SUSY-breaking 
matrix parameters m^,m^ and A,,, which govern the structure of the sneutrino squared- 
mass matrices. We have found a convenient parameterization of the sneutrino sector, where 
all relevant physical observables depend analytically on a pair of 3 x 3 mass matrices 
and given in eqs. (13.401) and (I3.39p . respectively. The elements of violate lepton 
number by two units, whereas elements of M^fj are lepton-number conserving parameters. 

Within this framework, we have analyzed the constraints arising from one-loop neutrino 
masses and mixings, from radiative flavor- changing charged lepton decays, and from the 
electron electric dipole moment (EDM). We discovered new and potentially significant 
contributions to radiative lepton decays — £^ + 7 due to the dependence of which 
modifies the MSSM value of M^^. We also observed that although the measurement 
places non-trivial constraints on the SUSY-breaking parameters, the electron EDMs do not 
yield any additional constraints (at one loop) on the seesaw-extended MSSM parameters. 
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All conclusions presented here are based on a complete numerical analysis of the processes 
described above ^ In all cases, we have also provided useful analytic approximations, which 
have served check of our numerical work. 

Sneutrino mixing phenomena takes on two different forms. The mixing of sneutrinos 
and antisneutrinos violates lepton number by two units, whereas sneutrino flavor mixing is 
a lepton-number conserving process. Both forms of mixing are in present in principle in the 
three-generation seesaw-extended MSSM. In this paper, we have generalized the sneutrino- 
antisneutrino mixing formalism, originally presented in a one-generation model [18], to 
the three-generation model. This sneutrino-antisneutrino mixing then acts back on the 
neutrino sector, and provides an important loop correction to the neutrino mass matrix. In 
this paper, we examined the possibility that starting from a diagonal neutrino mass matrix 
at tree-level, the nontrivial flavor structure of the neutrino mass matrix is generated entirely 
by the one-loop diagram that directly involves the sneutrino-antisneutrino transition. Our 
analysis shows that this is indeed possible, although in practice certain fine-tunings among 
SUSY breaking parameters in the leptonic sector seem to be unavoidable. 

Returning to the sneutrino sector, we have derived analytical expressions for both 
sneutrino-flavor oscillations and sneutrino-antisneutrino oscillations in eqs. (16.171) and 
f l6.18p . We determined that if the constraints analyzed above are combined with the as- 
sumption that sneutrinos can decay into two-body final states, then sneutrino-antisneutrino 
oscillations are not observable at colliders. This is consistent with a similar result of the 
one-generation model obtained in Ref. [18]. This conclusion is easily understood, by noting 
that the sneutrino-antisneutrino mass difference, Am,>, is proportional to the neutrino mass 
and is at most of the order of 1 keV. This is much smaller than the corresponding width 
of the sneutrino, Ty, of order 1 GeV or larger. The observability of sneutrino-antisneutrino 
oscillations at colliders requires that Amc, ~ Fj^. A sneutrino width of order 1 keV or less is 
possible only if there are no kinematically allowed two-body final states in sneutrino decay. 
In the seesaw-extended MSSM, this scenario is possible only if a charged slepton is the light- 
est supersymmetric particle, a possibility strongly disfavored by astrophysical data. Other 
possibilities exist if one introduces new degrees of freedom beyond the seesaw-extended 
MSSM, but this lies beyond the scope of this paper. 

Sneutrino flavor oscillations are more likely to be observable at colliders, since the 
mass splitting between sneutrinos of difference flavors can be of order 1 GeV or larger. 
We have derived simple approximate formulae for such oscillations and have estimated 
their magnitudes. Unfortunately, in the seesaw-extended MSSM, after imposing bounds on 
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Fortran-77 and Maple-10 numerical codes are available from the authors. 
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bounding sneutrino mixing angles determined from the analysis of radiative charged lepton 
decays, the resulting probabilities for sneutrino flavor oscillations are likely to be too small 
to be observed directly at colliders. 

At present, within the seesaw framework for neutrino masses, few handles exist for 
probing the physics at the seesaw scale. At most, one can hope to measure the MNS mixing 
angles, and determine neutrino mass differences (and with a little luck, the absolute scale of 
neutrino masses). In the seesaw-extended MSSM, some of the physics of the seesaw scale is 
imprinted on parameters that govern the properties of the light sneutrinos. With a precision 
program at future colliders for measuring sneutrino observables, there are new opportunities 
to explore the fundamental physics that is responsible for the origin of neutrino masses. 
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Appendix A Notation for fermion fields 



Fermion fields in quantum field theory can be described by employing either two-component 
or four-component fermion notation [56] . In models where lepton number is not conserved, 
two-component fermion notation is generally simpler and more efficient. In this appendix, 
we briefly discuss the relation between the two treatments. 

In Table [H the fermionic fields associated with the lepton and Higgs sectors of the 
seesaw-extended MSSM are listed. These fermion fields can be viewed either as two- 
component fermion fields or the left-handed projections of four-component fermion fields, 
with = 1(1 - 75)^ and 

^" = 0'^^, W=-^^C-\ (A.l) 

where = vl't-yO and C = —C^ is the charge conjugation matrix. 

For example, in four-component notation, given a four-component (anticommuting) 
Dirac spinor V£,, we define the following four- component spinors: 

vl = Pl^d , = Pl^d , = Pr^d , and = Pru'^ , (A.2) 

where Pl,r = |(1 respectively. The corresponding two-component (anticommuting) 

fields are given by the non-zero components of = Pl^d and z/£ = Pl^%- Consequently, 
we shall use the same symbols vl and for the corresponding two-component neutrino 
fields. However, one must be careful to note that in our notation 

vl = Cvf , ^ = -vie-' , (A.3) 



since, e.g., z/£ = PiCV^) = CIPrUd)^ . The same notation also applies to charged fermion 
fields. Our conventions for left and right-handed charged conjugated fields follow those 
of ref. [57]. Note that eq. flA.3l) implies that anticommuting fermion fields satisfy: 



^R^L = ^R^L , l^ll^R = l^Ll^R ■ (A.4) 

In the text, the effective Lagrangians for fermion mass and interaction terms are given in 
terms of two-component fermion fields. These terms can be easily translated into the four- 
component spinor notation . As a first example, the dimension-five operator that governs 
the standard seesaw mechanism [eq. (11.11) ] contains a product of two-component fermion 
fields, L^L^ . In terms of four-component spinors, this product is given by —{L'^)lC~'L^ = 
[R^y^L^ , where = {uf , i^) is now interpreted as a doublet of four-component fermion 
fields as described above and (-R'^)f = {i^r , ). 
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As a second example, we derive the four component version of eq. (13.11) in the one- 
generation model. One can redefine the phases of the neutrino fields such that mD and 
M are real and non-negative. The two-component spinor product z^l^£ + H.c. translates 
to the product of four- component spinors: —vJ^C'^vl + H.c. = UJiVl + WZun, which is the 
usual Dirac mass term. Similarly, the two-component spinor product z^£^'£ translates to the 
four-component spinor product —vf^C'^vl = T^i^l- Hence, if the Majorana mass term 
M 7^ in eq. (13. ip . one cannot identify the physical mass eigenstates as Dirac fermions. 
For example, the mass terms of the one-generation neutrino Lagrangian, which in terms of 
two-component fermion fields is given by — ^mass = ^d^l^l + j^^l^l + H.c, translates 
in four- component notation to 




(A.5) 



where we have used eq. (IA.4p to write the first line of eq. (lA.Sp in a symmetrical fashion 
and eq. flA.3p to obtain the final form above. 

The Takagi-diagonalization of the neutrino mass matrix yields two (self-conjugate) Ma- 
jorana fermion mass-eigenstates. This is accomplished by introducing a unitary matrix W, 

(A.6) 

) \ PlK J 

such that 

\U=\ - I , (A.7) 
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where m^^ ~ m^^/M and m,^,^ ~ M -|- m'jj/M. The resulting neutrino mass Lagrangian is: 

- J^mass = [m^.iyJC-'PLUe + m^^iyfC-'PLU^] + H.c. (A.8) 

We can define four-component self-conjugate Majorana fields by: 

ijM = Pl^i + PrCv^, , ?Af = ^^Pr - ^JC-^Pl , (A.9) 

= Pl^h + PRCvf , = KPr - J^fC-'PL . (A.IO) 
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Thus, eq. (1A.8I) reduces to the expected form: 



Appendix B A non-decoupling contribution to sneu- 

trino masses when ~ 0{M^) 

B.l Non-decoupling effects when ^ 

In Section [3l2| we noted below eq. (13.401) non-decoupling in the limit of — oo with 
||m^M~^|| fixed. The lepton-number conserving 3x3 squared-mass matrix of the light 
sneutrinos [eq. (13.391) ] can be written as: 

= ml + \Ml cos 2/5 + rrfj^M-^ml^M-^m^j^ + Oiy^M''^) + 0{v'^m%M-^) , (B.l) 

after expanding the quantity (1 + M~'^m%)^^ under the assumption that ||M~^m^|| < 1. 
Thus, we have a non-decoupling correction to the usual MSSM result of (9(m^M~^) as 
previously noted. 

To understand the origin of this non-decoupling phenomenon, we use eq. (15.41) which 
relates the original right-handed sneutrino with the light and heavy sneutrino states after 
block diagonalization of the sneutrino mass matrix. To formally integrate out the heavy 
sector and obtain the effective theory of the light sneutrinos, we must write: 

= Di - euniiM' + mlr'MYJY'LiHl , (B.2) 

before electroweak symmetry breaking, where we have used A^'^ = z/jj*. Note that when i7| 
is replaced by its vacuum expectation value V2/V^, we recover eq. (15.41) after using mo = 
V2Yy/\/2. In addition, we have used L^' ~ i>/ + 0{vM~^) and have worked consistently to 
leading order in vM~^. 

Consider the contribution of \dW / dN to the scalar potential, where W is given by 
eq. (El]). Then, 

^ = M^^N^ + e,,Yj<^H^Lf . (B.3) 

After squaring, and including the soft-SUSY-breaking term N*m\N (where is hermi- 
tian), we find: 

+ (z^) (^)* = ^^.--^"'^^"*^^^r£f^^* 

+ \tij{Y^,M)^^N^*H'lLf + H.c.l + (M^ + m^)^-^iV^*iV-^ . (B.4) 
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To obtain the relevant operator that survives in the low-energy effective theory, we insert 
eq. (1B.2I) for A^^ in eq. ( ]B.4I) . and then take the limit as ||M|| oo, In addition, we set 
i>h = 0. The end result is: 

e^nQ, [Y:yJ - Y:M{M^ + mlr'MYj] Li*LfHl*H^ . (B.5) 

Note that this is a dimension-4 (hard) SUSY-violating operator [58] which vanishes if 
m% = [as is the only SUSY-breaking source in eq. (1B.5I) ]. If m% < M^, one can 
expand (M^ + m^)"^ in eq. flB.5p . which yields: 

e,ne^Ay:M-'mlM-X^ + Oim%M-^)y^Li*Lf Hl*Hf . (B.6) 

We now replace vij V2. If ~ Oi^v"^), then the hard SUSY-breaking operator is of 

(9(f^M~^), which is the expected result. Such corrections are extremely small, assuming 
that V ^ and can be be dropped from the low-energy effective field theory of the 

light 0{v) degrees of freedom. On the other hand, if a; = ||m^||/||M^|| is held fixed to a 
finite positive value as M — oo, then the hard SUSY-breaking operator is of 0{x), which 
must be kept in the low-energy effective theory if x is not too small. 

In the latter case, we see the presence of a non-decoupling effect in the low-energy 
effective field theory of the 0{v) degrees of freedom as M — > oo. We identify this as a hard 
SUSY-breaking effect described by the dimension-4 operator given by eq. ( IB. 61) . Ultimately, 
this non-decoupling effect can be traced to the fact that although ul [I'l] and ul [z>^] are 
superpartners, it is not quite true that Ui [uh] and [uh] are superpartners. Explicitly 
[cf. eqs. (15. 3p and (15.41) ]. whereas 

u^,c^ul + M-'mluL, (B.7) 
to leading order in vM~^, we have: 

z>* ~z>^ + (M2 + m^)-iMm^i>i. (B.8) 

Clearly, with m% ^ 0, there is a slight discrepancy between Vh and the superpartner of Vh. 

If we replace li\ with its vacuum expectation value f 2/ a/2 in eq. ( IB. 51) and again make 
use of L{ ^ vf + C(t>M^^), we obtain a contribution to M\q: Then eq. (IB. 51) becomes: 

\m*^m^o - m*j^M{M^ + m\)-^ Mrn^^Y^ vi*vf , (B.9) 

which correctly reproduces the last two terms of given in eq. ( 13.391) . Of course, the 
non-seesaw MSSM result of MIq derives from the soft-SUSY-breaking term, L*m\Li, and 
the D-term contribution, |M| cos 2/5. As expected, in the M — 00 limit (with x ^ 0), the 
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low-energy effective theory reproduces the non-seesaw MSSM result. In this appendix, we 
have explained the origin of the non-decoupling correction to the non-seesaw MSSM result 
in the M — oo limit with x held fixed to a finite positive value. 

Finally, we address the question of the allowed size of the matrix parameter m^. Does 
it make sense to have x close to (9(1)? In ref. [38], it is shown that for values of x ~ 1, 
there is a very large negative shift in the mass of the lightest CP-even Higgs boson due 
to radiative corrections from the heavy neutrino/sneutrino sector of the seesaw-extended 
MSSM. If we demand that there should be no unusually large radiative correction to a 
physical observable generated as a result of 7^ 0, we can apply the results of ref. [38] 
for the radiatively-corrected physical Higgs masses to conclude that x < 0.1. Note that 
this upper bound is less severe than the bound of x < 0.01 given in eq. ( 14.18^ . The latter 
was obtained in Section 14.31 from the bounds on rare charged lepton radiative decay rates, 
which imply that the matrix should be close in form to a diagonal matrix. 

B.2 Naturalness constraints on the magnitude of 

It seems that phenomenological constraints allow for the possibility that ||?72,^|| is sig- 
nificantly larger than Olv"^), in which case the non-decoupling contribution to may be 
significant (perhaps as large as a few percent of the non-seesaw MSSM result). However, if 
one imposes the usual fine-tuning (or naturalness) requirements for the stability of the elec- 
troweak scale, one can show that cannot be significantly larger than (9(t'^). This can 

be verified by computing the one-loop correction to the if| self-energy. The computation 
in the supersymmetric limit is performed explicitly in Appendix E, section 7 of ref. [6] for 
the Wess-Zumino model. This computation is easily adapted to the present case of interest 
(in which the Higgs boson couples the the neutrino/sneutrino system). We then modify the 
supersymmetric computation in the case of the one-generation seesaw model by setting the 
boson (heavy sneutrino) squared-mass to + and the fermion (heavy neutrino) mass 
to M. [Here, we are dropping terms of (9(f^).] If 7^ (which softly breaks the super- 
symmetry), the quadratic divergence does not cancel exactly. The surviving contribution 
to the sqaured-mass term of if| is of the form 

m%\Y,\^I{M\m%)\H^W (B.IO) 

where T is a logarithmically divergent integral (that can be regularized by dimensional 
reduction [59]). 

We now add this one-loop result to the corresponding tree-level contribution to the 
scalar potential: 

iml^ + \fi\')\Hl\\ (B.ll) 
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In order to achieve successful electroweak symmetry breaking with v = 246 GeV, the 
complete coefficient multiplying must be of 0{v'^). By assumption, we take fi ~ 0{v) 
[cf. eq. f l2.12p ]. If ^ f^, the correct scale of electroweak symmetry breaking can 
be achieved only by an unnatural fine-tuning of the parameter m'jj^. Thus, naturalness 
requires that ~ v"^. We have not distinguished between 0{v'^) and O^M^^^y) in the 
above discussion. It is likely that there is a slight separation of scales with Msusy ^ 1 TeV. 
By imposing the naturalness condition on the dynamics of electroweak symmetry breaking 
(which ultimately is the motivation for TeV-scale supersymmetry in the ffist place), we 
conclude that the expected natural order of magnitude for is: 

llm^ll ~0(Msusy), (B.12) 

as indicated by eq. (12.151) . 

For completeness, we note that the same conclusion can be drawn by considering the 
one-loop effective scalar potential, V^^^cp). In particular, if we introduce a hard momentum 
cutoff A, one obtains a one-loop contribution of [60] 



^^'n</>) = ^ E StrM?^(0) + ^Str ^ Mf(0) 



In 



Mf(0) 1 
A2 2 



(B.13) 



where ((/)) are the contributing squared- mass matrices of particles whose masses originate 
from their couplings to the Higgs boson, with the vacuum expectation values replaced by 
the corresponding Higgs fields, 0, and 

Str {■ . ■ } = ^(-1)2^>(2J, + l)a {•■■}. (B.14) 



In eq. (1B.14|) . Ci counts the electric charge and color degrees of freedom of particle i (e.g., 
C = 2 for the gauge boson and C = 6 for a colored quark, since we count both particle 
and antiparticle). It is convenient to absorb the factor of 1/2 in the last term on the right 
hand side of eq. fIB.lSp . by defining /i such that: 

.n^P-i.ln^f^. (B.15, 

Using the results of eqs. (I3.12p . (13.691) and (13.701) . we focus on the contributions to the 
supertraces from the heavy neutrinos and sneutrinos. Indeed, 

^ Str {(f)) = 2 TTm% + 0{v^), (B.16) 

i 

although is field independent and thus contributes only to the vacuum energy. Here, 
we are interested in the implications of naturalness associated with electroweak symmetry 
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breaking (and not the cosmological constant). Thus we focus on the field-dependent part 
of the scalar potential that is quadratic in the Higgs fields. To do this, we simply replace 
niD with H2Y^. For simplicity, we shall examine the one generation seesaw model. In this 
case, we obtain the following scalar field-dependent squared-masses: 



ml ^ + 2\Y^\'\Hi\\ 



ml ^ M' + ml + \Y,\'\Hi\' 



M2 



M2 + 



m 



N 



(B.17) 
(B.18) 



Inserting these results into the last term on the right hand side of eq. fIB.lSp . and using 
eq. (IB.ISP to replace A with /i, we end up with the following terms in V^^^cp) that contribute 



to the coefficient of 



2 {(M^ + m%)^ + 2(2M2 + m%)\Y,\^\H^\^} In 
-2 {{M^ + AM^\Y,\^\Hi\^}\n 



r2 I _2 I IV |2| TJ212 ( 2M^+m% 



+ mi, + \Y,\'m 



M^ + 2\Y,\^\H^f 

7 , 



(B.19) 



where we have dropped terms of Expanding out the logarithms, the above 

expression reduces to 



2{{M^ + miy + 2{2M^ + ml)\Y,\^\Hi\^} \ In 



M'^ + m'^ 



N 



+ \YJ'\H. 



2| 7j2|2 



2M2 + 

^^2 



(M2 + m 



N) 



-2 {{M^ + 4M2|n|2|if2|2} (in ^ + 



2|y.p|g2T 



(B.20) 



If we keep only terms proportional to \H2\ \ we end up with: 



A%\'^\Hl\^ { 2MMn ( 1 + 



/r2 i„ / 1 , 



M2 



(B.21) 



One can check that the coefficient of is precisely m\1{M'^ ^m\), where X is the 



integral appearing in eq. fIB.lOp after DR subtraction [59] 



Appendix C Feynman rules 



We exhibit here the relevant Feynman rules for the calculation of 



?'7 presented in 



Section SSI These rules are based on four-component fermion notation (see Appendix A ) 
and employ the conventions of Ref. [7] for sfermion, chargino and neutralino masses and 
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mixing matrices. The neutrinos are (self-conjugate) Majorana fermions [cf. eq. flA.91) 
In the basis defined in Section [2] we obtain: 



J'-^MNS 



(C.l) 



—I 



g^Zi^Z]! - F/zi-'+^^^'z^* 



TjJI JD 



Lt-- 



Xi 



(C.2) 
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Lt-- 





A.I 



(C.3) 



(C.4) 



Appendix D Order of magnitude estimates for contri- 
butions to one-loop neutrino masses 

In this appendix, we estimate the order of magnitude of the one-loop contributions to the 
neutrino masses due to the graphs of fig. [2](a) and (b), and the corresponding graphs (not 
shown) in which the light sneutrinos [heavy neutrinos] in graph (a) [(b)] are replaced by 
heavy sneutrinos [light neutrinos]. 

In the case of graph (a), the dominant contribution involves the light sneutrino- 
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Figure 5: One- loop corrections to light neutrino masses. The x marks the location of 
the AL = 2 transition, (a) The loop consisting of light sneutrinos and gauginos. The x 
indicates the location of light sneutrino-antisneutrino mixing, and the solid dot indicates 
a factor of the gaugino Majorana mass in the numerator of the fermion-number-violating 
gaugino propagator, (b) The loop consisting of the neutral Higgs field and a heavy 
neutrino. The x indicates the lepton-number-violating heavy neutrino propagator, which 
is proportional to M6^^, and the solid dot indicates a mass insertion of the form 
The contributions of the corresponding graphs (not shown) in which the gauginos in (a) are 
replaced by the Higgsino -ff|, the light sneutrinos in (a) are replaced by heavy sneutrinos, 
and the heavy neutrinos in (b) are replaced by light neutrinos are all suppressed by an 
additional powers of 0{vM~^) as explained in the text. 

neutrino-gaugino interaction termQof eq. (15. 5p . We can estimate the leading contribution 
of this graph by replacing the internal lines by the interaction eigenstate fields that appear 
in eq. (15. 5p . as depicted in fig. [51 That is, we first replace the Sk with the z//, which must 
point away from both external vertices, as shown in fig. [5](a). The latter is possible only 
in the presence of light sneutrino-antisneutrino mixing, which is indicated by the x in 
fig. O^a). Using the expected magnitudes of the model parameters given by eqs. (12. lip and 
(I2.14p . the X in fig. [5](a) produces a factor Am|^ ~ 0{v^M~^). The neutralino line can be 
treated perturbatively. In the lowest order approximation, we take the neutralino to be a 
gaugino (either B or W^, with Majorana masses Mi and M2, respectively), and we treat 
the mixing of the gauginos with the neutral higgsino states {Hi and if|) as a perturbation. 
The corresponding gaugino propagators (with internal four-momentum q) shown in fig. [S](a) 
are fermion-number-violating propagators (indicated by the clashing arrows), and are given 
by iMk/{q^ — M|) for k = 1,2. We denote the presence of the gaugino mass [which is of 

^^Of the three hght sncutrino-neutrino-neutralino interactions of eq. (|5.5[) . the two sneutrino-neutrino- 
higgsino interaction terms are suppressed by a factor of ©(m^jAf^^) relative to the sneutrino-neutrino- 
gaugino interaction, and can be neglected. 
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0{v)] in the numerator by the sohd dot in fig. [5](a). Not including this exphcit factor of the 
gaugino mass, the loop in graph (a) then consists of two massive scalar propagators [with 
mass of 0(f )] and one fermion-number-violating propagator; hence the loop integral has a 
mass dimension of —2. Thus, the corresponding loop integral is of Combining the 

above results, the order of magnitude of the contribution of graph (a) is: 

which is indeed of order the tree-level neutrino mass multiplied by the product of the 
relevant vertex coupling constants and a typical loop factor of I/IGtt^ (denoted by Cl 
above). 

Suppose we replace the light sneutrinos of graph (a) with heavy sneutrinos. In this case, 
the effect of heavy sneutrino-antisneutrino mixing is Am?^ ~ 0{ml) ~ 0{vM). From 
eq. (15.51) . we see that there are potentially two contributions — one involving the gauginos 
and one involving the higgsino In the case of the gaugino loop graph, each vertex 
introduces a 0{vM~^) suppression. Thus, following the analysis above, we conclude that 
the order of magnitude of the heavy-sneutrino loop is suppressed by a factor of (9(t;^M~^) 
as compared with the light-sneutrino loop. In the case of the loop graph involving 
we note that there is no diagonal Majorana mass term for this higgsino field. Moreover, 
Hi does not couple to the external neutrinos, so we cannot use the off-diagonal Majorana 
mass term ^H\H2 for the fermion-number-violating neutralino propagator. Therefore, the 
heavy-sneutrino loop can be neglected. 

In the case of graph (b), the propagator of the heavy neutrino (with internal four- 
momentum q) is given by iM6^^/{q^ - M^), due to the presence of the lepton-number 
violating mass M (indicated by the x). Since the loop integral is dimensionless, it naively 
appears that the resulting loop integral should be of 0{M). However, an explicit compu- 
tation of the graph of fig. [2](b) demonstrates that the coefficient of the leading 0{M) term 
vanishes exactly after summing over the internal neutral Higgs and Goldstone states. The 
subleading term does not vanish and is of 0{v'^M^^), which is the magnitude of the light 
neutrino mass. This cancellation can be easily understood by noting that the two vertices 
of fig. [2](b) arise from interactions of eq. (15.51) that involve Thus we replace the neutral 
Higgs and Goldstone lines of fig. [2](b) by the if| field [cf. eq. (15. 6p ]. According to the 
interaction Lagrangian of eq. (15.51) . the if| field must point into both external vertices, as 
shown in fig. [5]^b). This requires a mass insertion on the if| line of the form (-ff|)^ + H.c. 
In fact, such a term exists in the MSSM Higgs potential [49] after shifting the neutral field 
H2 H2 + f 2/ -\/2, which results in a term of the form \'m\ sin^ /3(if|)^ + H.c. Thus, in the 
mass insertion approximation, graph (b) consists of the lepton-number- violating heavy neu- 
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trino propagator, two massive scalar field linecll and an insertion of 0{v'^). After extracting 
the factor of M from the numerator of the heavy neutrino propagator, the remaining loop 
integral now has a mass dimension of —2, which yields a result of 0{M~'^). Combining 
these result, the order of magnitude of the contribution of fig. [5](b) is given by: 

C'l^ ■M-v^ = C'J^^ , (D.2) 

which is again of order the tree-level neutrino mass multiplied by the product of the relevant 
vertex coupling constants and a typical loop factor (denoted above by C[). This result 
confirms our previous argument above. A careful evaluation of the leading behavior of the 
loop integral (in the limit of M ^ f ) then reproduces the result obtained in eq. (15. 8p . Note 
that the factor of sin (3 = f|/f2 that arises in the mass insertion on the if| line cancels 
out a similar factor of f| that appears in C[ oc Y^. 

If the heavy neutrinos in fig. [5]^b) are replaced by light neutrinos, the resulting contri- 
bution is suppressed by an additional factor of 0(t>^M~^) due to the suppression of the 
vlvf Hi interaction of eq. (15.51) . 
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